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ABSTRACT 

Studies  of  strength  development  at  polymer-polymer  interfaces  are  examined  and 
applications  to  welding  of  similar  and  dissimilar  polymers  are  considered.  The  fracture 
properties  of  the  weld,  namely,  fracture  stress,  a,  fracture  energy,  GIc,  fatigue  crack 
propagation  rate  da/dN,  and  microscopic  aspects  of  the  deformation  process  are 
determined  using  compact  tension,  wedge  cleavage,  and  double  cantilever  beam  healing 
experiments.  The  mechanical  properties  are  related  to  the  structure  of  the  interface  via 
microscopic  deformation  mechanisms  involving  disentanglement  and  bond  rupture.  The 
time  dependent  structure  of  the  welding  interface  is  determined  in  terms  of  the  molecular 
dynamics  of  the  polymer  chains,  the  chemical  compatibility,  and  the  fractal  nature  of 
diffuse  interfaces.  Several  experimental  methods  are  used  to  probe  the  weld  structure  and 
compare  with  theoretical  scaling  laws.  Results  are  given  for  symmetric  amorphous  welds, 
incompatible  and  compatible  asymmetric  amorphous  welds,  incompatible  semicrystalline 
and  polymer-metal  welds.  The  relevance  of  interface  healing  studies  to  thermal,  friction, 
solvent  and  ultrasonic  welds  is  discussed. 
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INTRODUCTION 

The  joining  or  welding  of  polymers  can  be  accomplished  by  several  techniques 
involving  thermal  welding,  vibrational  welding,  friction  welding,  solvent  welding,  surface 
chemical  modification,  ion  beam  surface  modification,  resonance  heating  and  other  more 
exotic  but  less  common  techniques.  In  this  Report  we  focus  on  polymer  welding  where  the 
molecules  near  the  surfaces  become  mobile  and  the  weld  strength  develops  by  a  combin¬ 
ation  of  surface  rearrangement,  wetting  and  diffusion.  The  molecular  mobility  can  be 
induced  by  heat  from  friction,  (ultrasonic)  vibration,  impact,  hot  plates,  solvents,  etc. 

In  recent  years,  we  have  considered  the  problem  of  strength  development  at  polymer- 
polymer  interfaces  in  terms  of  the  static  and  dynamic  properties  of  random-coil  chains  (1- 
4).  When  two  pieces  of  molten  polymer  are  brought  into  contact,  wetting  or  close 
molecular  contact  (van  der  Waals)  first  occurs  followed  by  interdiffusion  of  chain  segments 
back  and  forth  across  the  wetted  interface.  After  a  contact  time,  t,  we  inquire  as  to  the 
mechanical  energy,  G,  required  to  separate  the  two  pieces  as  a  function  of  time,  temper¬ 
ature,  T,  contact  pressure,  P,  and  molecular  weight,  M,  of  the  linear  random-coil  chains, 
i.e., 


G  =  W(t,  T,  P,  M)  (1) 

where  W  is  the  welding  function  to  be  determined.  Solutions  to  this  problem  have 
application  to  similar  situations  involving  polymer-polymer  interfaces.  For  example, 
processing  of  powder  and  pellet  resin,  internal  weld  lines  of  polymer  melts  during  extrusion 
and  injection  molding,  welding  of  surfaces,  lamination  of  composites,  coextrusion,  and 
autohesion  of  uncured  linear  elastomers. 

In  this  Report,  we  focus  on  problems  of  welding  and  autohesion  with  emphasis  on 
symmetric  interfaces  where  the  same  amorphous  polymer  exists  on  both  sides  of  the 
interface.  Consideration  is  also  given  to  the  important  case  of  incompatible  asymmetric 
amorphous  and  crystalline  interfaces  where  strength  develops  between  two  immiscible 
polymers.  We  begin  by  introducing  the  reptation  model  (5,  6)  for  diffusion  at  a  symmetric 
interface  and  determine  a  molecular  description  of  the  interface  structure  as  a  function  of 
the  variables,  t,  T,  P,  and  M.  We  then  discuss  models  to  relate  the  interface  structure  to 
the  mechanical  properties  via  a  set  of  microscopic  deformation  mechanisms.  The  latter 
involve  disentanglement  by  chain  pullout,  bond  rupture  and  connectivity  between  chains. 
The  predictions  for  welding,  tack  and  green  strength  are  compared  with  experiments  and 
theories  (7-11)  of  other  investigators. 
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THEORY 


Molecular  Dynamics  of  Random-Coil  Chains 

The  motion  of  individual  chains  in  amorphous  bulk  materials  or  concentrated  solutions 
of  linear  random-coil  (Gaussian)  chains  has  been  modelled  by  the  reptation  theory  of  de 
Gennes  (5)  and  Edwards  (6).  In  this  model  the  chain  is  confined  to  a  tube  having  a  similar 
shape  as  the  random-coil  configuration  of  the  chain.  The  tube  represents  topological 
constraints  to  la:  'ral  motion  of  monomers  imposed  by  the  neighboring  chains  via  entangle¬ 
ments  and  predc  ninantly  restricts  motion  of  the  chain  to  that  along  the  curvilinear  length 
of  the  tube.  Our  use  of  this  model  is  illustrated  in  Fig.  1.  At  t  =  0,  the  chain  (solid  line)  is 
in  its  initial  tube  indicated  by  the  dashed  lines.  The  chain  exhibits  Brownian  motion  back 
and  forth  in  the  tube  and  because  the  chain  ends  are  free  to  move  in  any  direction  away 
from  the  tube,  the  memory  of  the  initial  tube  position  in  space  is  gradually  lost  as  shown  at 
times  t2  and  ti.  However,  the  chain  being  in  a  dense  system  is  always  in  a  new  tube  but 
only  parts  of  tne  chain,  usually  towards  the  center,  retain  memory  of  the  initial  tube  for  the 
relaxation  time,  Tr.  The  length  of  chain,  l(t),  which  "escapes”  from  the  initial  tube  is  also  a 
random-coil  chain  obeying  Gaussian  statistics  and  we  call  it  the  minor  chain  (2).  As  the 
minor  chain  increases  in  length  we  can  represent  its  static  dimensions  by  a  "most  probable" 
spherical  envelope  as  shown  in  Fig.  1. 


0<  1,  <  12<  Tr 


Fig.  1.  The  reptation  model  for  a  random-coil  chain  in  an  entangled  melt  shows  the  chain 
relaxing  from  its  original  configuration  or  tube  at  times  Octj  <t2^TV  where  Tr  is  the  tube 
disengagement  or  relaxation  time.  The  minor  chain  is  the  portion  of  the  chain  which  has 
relaxed  from  its  original  configuration  and  is  shown  contained  in  its  random-coil  most 
probable  spherical  envelope.  The  dotted  tube  represents  the  topological  constraints  to 
motion  imposed  by  the  neighboring  chains  in  the  melt  (Kim  and  Wool). 


Finally  at  t  =  Tr,  the  reptation  time,  the  chain  has  escaped  or  forgotten  its  original 
configuration.  The  following  relations  summarize  the  description  of  the  motion  of  a  linear 
polymer  chain  in  an  entangled  melt: 


Tlie  mean  square  escape  length  of  the  minor  chains,  <12>: 


<-l2>  =  16Di0r. 

The  curvilinear  one-dimensional  diffusion  coefficient,  Dj: 


(2) 
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Dj  «  M'1.  (3) 

The  mean  square  monomer  displacement  of  the  minor  chains: 

<X2>  *  <12>  V2.  (4) 

The  mean  square  center-of-mass  diffusion  distance,  <  X2Cm  > : 

<^2cm>  ~  2Dt.  (5) 

The  center-of-mass  self-diffusion  coefficient,  D: 

D  -  M~2.  (6) 

The  reptation  time,  Tr: 

Tr  -  M3  (7) 

The  fraction  of  chain  remaining  in  the  tube,  F(t): 

F(t)  =  1  -  4/JI3/2(t/Tr)1/2  (8) 

where  t  <  Tr.  For  t  >  Tr,  F(t)  can  be  well  approximated  by  (5,  12), 

F(t)  ~  8/jr2e“t/^'r.  (9) 


Equations  8  and  9  form  the  basis  for  describing  the  viscoelastic  shear  relaxation 
modulus  of  the  Doi-Edwards(12)  constitutive  theory  of  entangled  linear  polymer  melts. 
Considerable  experimental  support  exists  for  the  molecular  weight  dependence  of  the  self¬ 
diffusion  coefficient,  Eq.  5,  and  has  been  reviewed  by  Tirrell  (13).  Support  for  the  minor 
chain  model  shown  in  Fig.  1  has  recently  been  obtained  by  Lee  and  Wool  and  by  Walczak 
and  Wool  (14).  Using  FTIR  techniques,  they  showed  that  the  orientation  relaxation  of 
step-strained  centrally  deuteriated  polystyrene  chains  first  occurred  by  relaxation  of  the 
protonated  chain  ends  while  the  center  of  the  chain  retained  its  orientation.  The  chains 
were  diluted  in  a  higher  molecular  weight  matrix.  The  centrally  deuteriated  part  of  the 
triblock  chain  then  relaxed  at  a  predicted  time  and  rate  consistent  with  the  reptation 
theory.  The  homopolymer  chains  (without  the  higher  molecular  weight  matrix)  do  not 
show  this  effect  which  we  believe  is  due  to  the  percolation  relaxation  mechanism  recently 
proposed  by  us.  However,  the  reptation  theory  fails  to  accurately  predict  the  molecular 
weight  dependence  of  the  zero-shear  viscosity  with  a  3.4  power  law  exponent.  The 
relaxation  time,  Tr  «  M  ,  is  also  important  in  determining  healing  or  welding  times  at 
polymer  interfaces.  In  the  next  section  we  use  the  minor  chain  model  to  develop  a  time 
and  molecular  weight  dependent  molecular  description  of  a  polymer  interface.  Effects  of 
temperature  and  pressure  will  also  be  discussed. 

Molecular  Description  of  a  Polymer-Polymer  Interface 

The  most  important  interface  to  consider  is  that  formed  by  a  symmetric  A/A  polymer 
pair  consisting  of  linear  monodisperse  molecular  weight  polymers  whose  chain  ends  are 
uniformly  distriouied  in  the  bulk  melt  material.  Polydisperse  symmetric  A/A  pairs  and 
some  asymmetric  A/B  pairs  can  also  be  investigated  using  methods  presented  in  this 
section.  Assuming  instantaneous  wetting  (time-dependent  wetting  is  considered  later)  of 
the  two  surfaces  we  proceed  to  analyze  the  molecular  aspects  of  the  interface  arising  from 
interdiffusion.  This  problem  is  simplified  further  using  the  minor  chain  model  as  shown  in 
Fig.  2. 


(o) 

INTERFACE 
I  .  *  CNH  SEGMENT 


(b) 


Fig.  2.  The  interdiffusion  process  at  a  polymer-polymer  interface  in  terms  of  the  behavior 
of  (a)  the  minor  chains  and  (b)  the  minor  chain  spherical  envelopes.  Only  one  side  of  the 
interface  is  shown  for  convenience.  Instantaneous  molecular  contact  or  wetting  is  assumed 
at  t=0  (Kim  and  Wool). 

At  t  =  0,  molecular  contact  is  achieved  at  the  interface  (only  the  molecules  on  one  side 
of  the  interface  are  shown  for  clarity).  The  chains  begin  to  reptate  via  their  chain  ends 
which  are  shown  initially  as  dots  randomly  distributed  in  space.  As  time  proceeds,  only  that 
part  of  the  chain  which  has  escaped  its  original  configuration,  i.e.,  the  minor  chain,  can 
contribute  to  interdiffusion  across  the  interface.  Furthermore,  we  do  not  have  to  consider 
the  exact  location  of  each  segment  of  the  minor  chain  but  only  the  spherical  envelope  to 
which  it  belongs,  as  shown  in  Fig.  2.  Thus,  as  time  approaches  the  reptation  time,  Tr,  inter¬ 
penetration  and  reentanglement  of  chains  are  achieved  and  the  molecular  properties  of  the 
virgin  state  are  obtained  at  the  interface  plane.  From  this  minor  chain  model,  many 
molecular  aspects  of  the  interdiffusion  process  can  be  determined,  as  presented  below. 
Solutions  are  given  for  t  <  Tr  (healing)  and  t  >  Tr  (virgin  state).  Details  of  the 
calculations  are  given  in  references  2,  3,  and  15. 


Number  of  Chains  Intersecting  the  Interface,  n(t) 

The  number  of  random-coil  chains  intersecting  unit  area  of  the  interface  as  a  function 
of  contact  time,  t,  and  molecular  weight,  M,  behaves  as, 

n(t)  *  t1/4  M'5/^;  (t  <  Tr)  (10) 

noo  «  M"1/2.  (t  >  Tr)  (11) 

A  chain  is  considered  to  intersect  the  interface  when  one  of  its  ends  has  crossed  the 
interface  contact  plane  defined  at  t  =  0.  n(t)  also  scales  as  X/M  where  X  is  the  average 
monomer  segment  interpenetration  distance  defined  by  Eq.  4.  Eq.  10  is  derived  on  the 
basis  that  all  chains  whose  ends  are  within  a  distance  of  X  from  the  interface  will  have 
intersected  the  interface  at  time  t,  and  that  the  number  of  chadns  per  unit  volume  varies  as 
M  .  The  time-dependence  is  a  typically  t1/4  rather  than  t1/2  as  expected  for  usual  atomic 
diffusion  processes. 

Note  that  the  staticpr  virgin  state  solution  at  t  >  Tr  is  obtained  by  substituting  for  the 
reptation  time,  Tr  ~  in  the  dynamic  solution.  The  static  solutions  are  independent  of 
the  molecular  dynamics  model  and  can  be  derived  from  other  considerations.  Eq.  11  was 
derived  previously  by  Wool  and  Rockhill  (17).  Thus,  the  dynamic  solutions  give  the  static 
solutions  simply  by  substituting  for  the  molecular  weight  dependence  of  the  relaxation 
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time.  This  is  an  important  feature  in  evaluating  molecular  aspects  of  strength 
development. 

When  a  molecular  property  is  suggested  as  controlling  the  time  dependence  of  welding, 
that  property  can  be  used  to  simultaneously  predict  both  the  molecular  weight  dependence 
of  welding  and  the  molecular  weight  dependence  of  fracture  in  the  virgin  state.  This  three 
way  concurrence  is  useful  in  evaluating  potential  solutions  to  welding  problems.  For 
example,  it  is  known  from  the  experiments  of  Kausch,  et  al.  (18-20)  and  Wool,  et  al.  (1,  21, 
22)  that  the  fracture  stress  during  crack  healing  and  welding  increases  as  a  ~  tr/  .  It  could 
be  hypothesized  upon  inspection  of  Eq.  10  that  the  number  of  chains  intersecting  the 
interface,  n(t)  ~  t1/  ,  is  responsible  for  the  strength  development.  However,  this 
hypothesis  can  be  rejected  by  examining  the  virgin  state  property,  Eq.  1 1,  since  the  fracture 
stress  of  polymers  does  not  decrease  with  increasing  molecular  weight  as  Ooo  ~  M*  V2 
Further  experimental  study  would  show  that  the  healing  rate  does  not  depend  on  M"5'  . 
Conversely,  molecular  properties  used  to  describe  the  virgin  state  strength  can  also  be 
evaluated  by  their  ability  to  describe  the  time  and  molecular  weight  dependence  of  healing. 
This  argument  applies  when  the  same  microscopic  deformation  mechanisms  are  operative 
during  welding  and  fracture  of  the  virgin  state. 


Average  Interpenetration  Contour  Length,  l(t) 

The  average  contour  length  of  chain  segments  or  cilia  which  have  diffused  across  the 
interface  is  obtained  from  the  minor  chain  model  as, 


l(t)  *  t1/2  M'1/2  (12) 

loo  *  M  (13) 

This  property  has  the  same  scaling  relation  as  the  minor  chain  length,  <  1  > .  It  will  be 
shown  later  that  l(t)  is  important  in  determining  mechanical  properties  where 
disentanglement  and  chain  pullout  are  dominant  microscopic  deformation  mechanisms 


Number  of  Bridges  Crossing  Interface,  p(t) 

Each  time  a  piece  of  interdiffusing  chain  crosses  back  and  forth  through  the  interface, 
it  creates  a  molecular  bridge.  Intuitively,  this  may  be  visualized  as  a  "sewing  up"  of  the 
interface.  p(t)  is  determined  by  p(t)  =  n(t)Vl(t)  where  Vl(t)  is  proportional  to  the  number 
of  bridges  created  by  each  minor  chain  (a  random  walk  of  N  steps  originating  at  a  plane 
will  cross  the  plane  VN  times).  Thus  we  obtain, 

p(t)  =  t^/2  M'2/2  (14) 

Poo  *  M°.  (15) 

This  result  was  first  obtained  by  de  Gennes  (9)  using  a  scaling  law  and  is  equivalent  to 
the  "crossing  density"  determined  by  Prager  and  Tirrell  (7)  from  a  more  detailed  calcula¬ 
tion  of  the  mterdiffusion  process.  Note  that  the  virgin  state  number  of  bridges  is 
independent  of  molecular  weight.  Thus,  bridges  alone  cannot  be  used  to  describe  strength 
development  during  welding  in  cases  where  the  virgin  strength  exhibits  a  molecular  weight 
dependence  over  some  range  of  molecular  weight  as  observed  for  most  polymers. 

Modifications  to  the  bridge  models  have  been  suggested  by  Adolf  Tirrell  and  Prager  (7) 
and  later  by  Mikos  and  Papas  (29)  which  involve  a  minimum  interpenetration  depth,  Xe,  to 
develop  entanglements.  These  approaches  introduced  a  molecular  weight  dependence  for 
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the  virgin  state,  resolving  the  problem  posed  by  Eq.  15  and  gave  a  more  complex  expression 
foi  the  time  dependence  of  healing.  Extrapolation  of  long  time  healing  data  to  zero 
strength  is  predicted  by  this  modification  to  intersect  th°  time  axis  at  a  finite  time,  t^, 
required  to  diffuse  to  the  minimum  distance  Xe.  This  effect  has  not  been  observed  in  many 
welding  studies  conducted  by  various  groups. 

Polydispersity  effects  were  also  considered  by  Adolf,  Tirrell  and  Prager  (7).  Assuming 
independent  reptating  chain  diffusion,  the  contribution  to  the  concentration  profile  from 
individual  chains  can  be  readily  determined  for  any  molecular  weight  distribution.  At  short 
times,  the  contribution  from  the  faster  diffusing  shorter  chains  dominates  the  weld  strength 
followed  later  by  the  slower  diffusing  longer  chains.  This  can  result  in  highly  nonlinear 
welding  behavior  as  demonstrated  by  McGarel  and  Wool  (54)  for  welding  of  polystyrene. 
Bridge  concepts  are  also  useful  in  the  development  of  entanglement  models  and  for 
estimating  the  number  of  broken  bonds  in  brittle  fracture. 


Average  Monomer  Interpenetration  Depth,  X(t) 

The  average  monomer  interpenetration  depth  of  segments  on  the  interdiffusing  chains, 
X(t),  is  obtained  from  an  integration  over  the  Gaussian  segment  density  profiles  of  the 
interdiffused  minor  chain  population,  (2)  as  shown  in  Fig.  3,  such  that, 

X(t)  ~  t1/4  M'1/4;  (16) 

Xoo  -  mV2  «  R/V6  (17) 

At  t=Tr,  the  average  monomer  interpenetration  distance  is  equal  to  the  radius  of 
gyration,  or  X =R/V6.  The  t*/4  dependence  in  Eq.  16  is  unique  to  the  correlated  motion 
of  the  chains  for  times  less  than  the  reptation  time.  This  result  can  be  understood  by 


Fig.  3.  Segment  density,  p,  of  minor  chains  and  the  average  interpenetration  distance,  x 
(Kim  and  Wool). 

considering  the  Brownian  motion  and  displacement,  s,  of  an  individual  monomer  along  the 
contour  of  the  tube  as  s  ~  'for  t  <  Tr.  However,  the  curvilinear  motion  is  related  to 
translational  motion  via  X  ~  s1/^  because  the  chain  is  a  random-coil.  For  times  greater 
than  Tr,  when  the  original  tube  has  evaporated  the  motion  of  the  monomer  becomes 
uncorrelated  since  s  is  greater  than  the  chain  length  and  its  displacement  becomes  similar 
to  the  center-of-mass  displacement,  Xcm, 


Xcm  »  tV2M-l 


(18) 
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The  crossover  from  t*/4  to  behavior  could  be  examined  by  experimental 
techniques  such  as  Secondary  Ion  Mass  Spectrometry  (SIMS),  Sputtered  Neutral  Atom 
Mass  spectrometry  (SNMS),  (23)  and  Neutron  Specular  Reflectivity  (NSR),  which  probe 
the  concentration  profile  at  a  polymer  interface  for  t  <  Tj-  and  distances  X  <  Rp.  Center 
of  mass  motion  for  t  >  Tr  and  X  »  Rp  is  well  established  from  self-diffusion  studies  and 
found  to  be  in  agreement  witn  Eq.  1£  For  chains  in  the  bulk,  the  monomer  motion  is 
predicted  to  proceed  as  t1/4  for  t  <  Tr  while  the  center-of-mass  proceeds  as  t1/2  at  all 
times.  Doi  and  Edwards  (12)  showed  that  the  center-of-mass  motion  is  uncorrelated  at  all 
times  by  virtue  of  the  random  motion  of  the  chain  ends. 

Eq.  18  cannot  be  used  to  describe  center-of-mass  motion  for  chains  whose  random-coil 
configurations  a;e  altered  by  the  interface  as  shown  in  Fig.  4.  The  chains  initially  have  a 
non-Gaussian  configuration  due  to  the  reflecting  boundary  condition  at  the  surface.  As  the 
minor  chains  interdiffuse,  the  configurations  of  the  affected  chains  in  the  surface  layer  relax 
to  Gaussian  configurations.  In  Fig.  4,  the  average  shape  and  size  of  a  chain  at  the  interface 
(solid  vertical  line)  before  and  after  relaxation  are  shown.  The  distance,  77,  from  the 
location  of  the  highest  segment  density  of  the  chain  to  the  interface  remains  constant 
during  the  relaxation.  Since  the  center-of-mass  of  each  chain  is  close  to  the  peak  of  the 
segment  density  profile,  it  can  be  seen  that  v*ry  little  motion  of  the  center-of-mass  occurs 
for  t  <  Tr.  For  a  chain  in  the  bulk,  the  center-of-mass  is  displaced  by  a  distance  approxi¬ 
mately  equal  to  its  end-to-end  vector,  R,  during  iie  same  relaxation  time.  Since  R  >  tj,  the 
"apparent  center-of-mass  self-diffusion  coefficient"  of  chains  affected  by  the  interface 
would  be  smaller  than  that  of  chains  further  away  from  the  interface  in  the  bulk.  This 
argument  underlies  the  qualification  for  the  center-of-mass  scaling  law  in  Table  1. 

The  nonecpiilibrium  configurations  of  the  chains  near  the  interface  plays  an  important 
role  in  determining  entanglement  density  and  the  structure  of  asymmetric  incompatible 
amorphous  interfaces.  The  reduced  entropy  of  the  confined  chain  configurations  shown  in 
Fig.  4  acts  as  a  driving  force  to  achieve  limited  diffusion  in  the  interface  of  immiscible 
polymers.  The  increase  in  entropy  is  counterbalanced  by  a  positive  enthalpy  of  mixing  and 
results  in  an  equilibrium  thickness  predicted  by  Helfand,  et  al.  (24)  as  doo  ~  Chi'1'72  where 
Chi  is  the  Flory-Huggins  interaction  parameter  (25). 

to) 


Fig.  4.  (a)  Average  shape  and  size  of  a  chain  at  the  interface  before  and  after  healing,  (b) 
Segment  density,  p,  of  the  chain  in  (a)  before  and  after  healing  (Kim  and  Wool). 


Number  cf  Monomers  Crossing  Interface,  L0(t) 


The  number  of  monomers  which  have  crossed  from  one  side  of  the  mterface  to  the 
other  is  equivalent  to  the  total  interpenetration  contour  length,  L0(t)  =  n(t)l(t).  Summing 
over  the  segment  density  profiles  of  the  interdiffused  minor  chains,  we  obtain,  (15) 

Lc(t)  -  t3/4  M'7/4;  (19) 

LO00  -  M1/2.  (20) 

The  unusual  t3/4  dependence  :s  again  due  to  correlated  motion  effects  and  reverts  to  a 
normal  t1/^  diffusion-dependence  for  times  greater  than  Tr.  This  result  has  been 
supported  by  experimental  techniques  discussed  in  reference  23  where  L0  is  determined 
from  the  integral  of  the  concentration  profile. 


Concentration  Profile 

Integrating  over  the  population  of  interdiffused  minor  chain  Gaussian  segment  density 
profiles  shown  in  Figs.  2  and  3,  the  normalized  monomer  concentration  profile  <X>(x,t),  for 
the  symmetric  polymer-polymer  interface  was  obtained  by  Zhang  and  Wool  (15)  for  t  <  Tr 

as, 


<b(x,t)  =  [l(t)/L  +  x2/bL]  erfc[x/(bV2n)]  (21) 

-  V(2n/jr)  (x/L)  exp  [-x2/(2nb2)] 

where  b  is  the  bond  length,  n  is  the  number  of  monomers  in  the  minor  chain  of  length  l(t), 
x  is  the  diffusion  depth  such  that  x  =  0  at  the  interface  and  L  is  the  contour  length  of  the 
entire  chain.  The  time  dependence  of  the  concentration  profile  in  Eq.  21  derives  from  the 
time-dependence  of  the  minor  chains  via  I(t)  =  V(16  Dj  0r). 

Computed  concentration  profiles  are  shown  in  Fig.  5  at  fractional  reptation  times  and 
diffusion  depths  normalized  with  respect  to  the  end-to-end  vector  of  the  Gaussian  chains 
(151.  Fig.  5  clearly  shows  that  the  concentration  at  the  initially  joined  interface  at  x  =  0  is 
discontinuous  when  the  welding  time  is  shorter  than  the  reptation  time.  For  normal 


Fig.  5.  Concentration  profile  of  a  polymer-polymer  interface  at  timer  t  less  than  the 
reptation  time,  Tr.  The  interdiffusion  depth,  x,  is  normalized  with  respect  to  the  radius  of 
gyration,  Rg  (Zhang  and  Wool). 


Fickian  dinusion,  the  atomic  or  monomer  concentration  profile  at  the  origin  immediately 
attains  a  value  of  0.5,  (assuming  an  initially  sharp  interface  at  t  =  0)  since  the  Brownian 
motion  of  the  atoms  is  isotropic.  For  polymers  howevei,  fhe  majoritv  of  monomers  in 
contact  with  the  initial  plane  exhibit  anisotropic  motion  due  to  their  connectivity.  1  hey  are 
not  free  to  interdiffuse  on  the  time  scale  of  monomer  motion  but  await  the  diffusion  of  the 
much  smaller  number  of  minor  chain  ends.  Thus,  <I>(x,t)  contains  a  discontinuity  at  x  =  0 
which  is  an  interesting  peculiarity  of  the  reptation  model  and  persists  for  at  least  one 
reptation  time.  The  discontinuity  or  gap  may  be  difficult  to  observe  experimentally 
because  the  short  range  segmental  mot.on  may  obscure  the  gap  rear  x  =  0  (15).  However, 
analysis  of  the  entire  profile  shape,  particularly  at  x  »  0,  should  give  crit’^al  insight  into  the 
reptation  model.  This  work  is  in  progress  using  Neutron  Reflection. 

For  t  >  Tr,  the  concentration  profile  returns  to  the  classical  Fickian  profile, 

<F(x)  =  erfc(x/L(j),  (22) 

where  is  the  Einstein  diffusion  length  given  by  =  2\^(Dt).  The  non-Fickian 
characteristics  of  the  concentration  profile  due  to  correlated  motion  effects  in  the  reptation 
model  have  been  discussed  by  de  Gennes,  (9)  Prager  and  Tirrell,  (28)  and  by  Zhang  and 
Wool  (  i5).  Eq.  21  has  also  been  used  to  calculate  the  molecular  properties  of  the  interface 
~nd  these  results  are  summarized  in  Table  1  (15). 

Scaling  Laws  for  a  Polymer-1  olymer  Interface 

The  molecular  properties  described  above  and  others  listed  in  Table  1  have  a 
convenient  common  scaling  law  which  relates  the  dynamic  properties,  H(t),  to  the  static 
properties,  Hoc,  via  the  reduced  time,  t/Tr,  by  (4) 


H(t)  =  Hoo  (t/Tr)r/4; 

(23) 

Hoc  -  M(3r’s)/4, 

(24) 

where  r,  s  =  1,2,3, . 

Table  1  summarizes  the  average  dynamic  and  static  molecular  properties  of  a  polymer- 
polymer  inteiface  and  lists  values  of  i  and  s  for  individual  properties.  This  table  can  be 
used  to  determine  other  properties  and  evaluate  suggestions  and  theories  relating  interface 
structure  to  mechanical  properties.  Several  of  these  properties  have  also  been  derived  by 
de  Gennes  (9,  26,  27)  and  by  Prager  and  Tirrell  (7,  28,  29). 


Table  1.  Molecular  Aspects  of  Interdiffusion  at  a  Polymer-Polymei  Interface. 


Molecular  Aspect 

Symbol 

Dynamic  Relation 
H(f) 

Static  Relation 

H. 

r 

$ 

General  property 

Hit ) 

t'l‘M~‘l* 

r 

s 

No.  of  chains 

nit) 

/’"Ar5" 

1 

5 

No.  of  bndge 

Pit) 

2 

6 

Ave.  monomer  depth 

X(t) 

Mxrl 

1 

1 

Total  monomer  depth 

Xo(t) 

M° 

2 

6 

Ave  contour  lengtt, 

fit) 

t'nM'r‘ 

M 

2 

2 

Total  contour  length,  no.  of  monomers,  N 

Ut) 

M'12 

3 

7 

Avr  bndge  length 

Ut) 

t  'I'M''1* 

M'13 

1 

1 

‘  Center  of  mass  depth 

M'V 

2 

5 

Diffusion  frot.  'ength 

N, 

|0 

2 

6 

*  TtW*  aquation  appka«  to  chain#  In  tha  bulk  and  doe*  not  apply  to  chakra  whoaa  cantar  of  mail  It  wtthkr  a  radio*  of  gyration  «f  tha  turtaca 
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If  a  mechanical  property  such  as  the  fracture  energy,  G,  is  controlled  by  one  of  the 
properties,  H(t),  then  the  time-depenence  of  welding  is  given  by  G  ~  tr' 4  and  the 
molecular  weight  dependence  by  M's/  .  Considering  the  complexity  of  the  fracture 
mechanics  and  deformation  mechanisms,  it  may  be  an  overly  simplified  argument  to 
suggest  that  a  single  molecular  property  defined  typically  over  a  few  hundred  angstrom 
thickness  of  the  interface  is  responsible  for  the  energy  dissipation  and  behavior  of  a  much 
larger  deformation  zone  at  the  crack  tip.  However,  the  scaling  laws  presented  in  Table  1 
offer  a  convenient  framework  for  evaluating  molecular  models  in  terms  of  the  static  and 
dynamic  properties  of  the  polymer  chains. 

FRACTAL  NATURE  OF  A  DIFFUSION  FRONT 


Atomic  Diffusion  Front 

While  the  concentration  profile  due  to  diffusion  varies  smoothly  with  depth  as  shown  in 
Fig.  5,  the  diffusion  field  when  viewed  in  two  or  three  dimensions  has  a  veiy  rough  nature. 
Consider  the  two-dimensional  lattice  diffusion  of  A-atoms  (Fig.  6)  into  B-atoms  (not 
shown)  as  discussed  by  Sapoval,  et  al  (30-32).  The  concentration  profile  of  A-atoms  is 
given  by  the  Fickian  result,  <F(x)  =  erfc(x/Ld),  where  L(j  =  2(Dt)^/^  is  the  Einstein 
diffusion  length  at  time,  t,  with  diffusion  coefficient,  D.  The  diffusion  "front"  (heavy  line  in 
Fig.  6)  is  defined  by  those  A  atoms  which  are  connected  to  the  diffusion  source  at  X  =  0 
via  other  A  atoms  and  have  an  empty  first  or  second  B  neighbor  which  is  itself  connected  to 
the  source  of  B  atoms.  The  interface  with  a  diffusion  front  is  seen  to  consist  of  solid 
regions  with  lakes  or  holes  leading  down  to  the  "seashore"  and  unconnected  islands  of  A 
atoms  in  the  B  "sea".  The  presence  of  holes  and  the  connectednet  oetween  atoms  in  the 
interface  is  important  in  determining  electrical  and  mechanical  pi  jperties  of  the  interface. 


A  Atom 


Fig.  6.  Schematic  of  a  two-dimensional  diffusion  field  at  a  given  time.  The  heavy  line 
represents  the  "Diffusion  Front"  determined  by  the  leading  edge  of  A-atoms  connected  to 
their  diffusion  source  at  x  =  0  and  having  at  least  one  B-neighbor  connected  to  its  source. 
The  position  of  the  center  of  the  front  occurs  at  Xf  and  the  width  of  the  front  is  2^f  (From 
Sapoval,  et  al.,  Reference  30). 

Sapoval,  et  al.  (30-32),  determined  the  frontier  to  have  fractal  properties  involving  self¬ 
similarity  (scale  invariance)  and  incompleteness  of  space  filling  (holes).  The  fractal 
geometry  of  the  frontier  is  determined  by  the  mass  to  radius  relation,  (33) 
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M  «  RDf,  (25) 

where  Df  is  the  fractal  dimension  of  the  diffusion  front.  The  fractal  geometry  exists  only 
over  the  width  of  of  the  front  as  shown  in  Fig.  6.  In  two  dimensions.  Dr  =  1.76.  The 
number  of  particles  per  unit  width  of  the  front,  Nf,  increases  with  the  diffusion  length 
according  to 


Nf  =  0.96  Lda,  (26) 

where  the  exponent  a  is  given  by, 

a  =  (Df-l)v/(l+v),  (27) 

in  which  v  =  4/3  is  the  critical  exponent  for  the  percolation  coherence  length  in  two 
dimensions,  resulting  in  a  =  0.43.  The  latter  result  is  also  supported  by  computer 
simulation  studies  of  Wool  and  Long  (34).  The  mean  position  of  the  frontier  can  be  shown 
to  occur  at  a  concentration  corresponding  to  the  percolation  threshold,  pc  (30).  In  two 
dimensions,  pf  =  0.59,  such  that  particles  in  the  profile  at  3><pc  are  not  accessible  or 
connected  while  particles  in  the  gradient  at  <h>pc  are  connected  to  the  source  at  x  =  0. 

The  three-dimensional  analysis  of  diffusion  fronts  is  more  complex  (31). 

Eq.  26  indicates  how  the  interface  "roughness"  increases  with  diffusion  length  such  that 
Nf  ~  L(iu'  .  Thus,  even  though  the  average  concentration  profile  varies  smoothly,  the 
diffusion  front  can  have  considerable  lateral  roughness  as  shown  in  Fig.  6.  Experimental 
examples  of  this  behavior  are  described  in  reference  34  based  on  studies  of  electrochemical 
reduction  of  metal  ions  diffusing  in  polymers  by  S.  Mazur  (35).  Fig.  7  shows  a  polymer 
metal  interface  produced  by  electrochemically  depositing  silver  atoms  in  a  polyimide  film. 


MONOMER  INTERFACE 


Fig.  7.  Computer  simulation  of 
(a)  monomer  interface,  (b)  polymer 
interface,  and  (c)  the  connected 
region  of  a  silver/polyimide  inter¬ 
face  (Wool  and  Long). 


POLYMER  INTERFACE 


POLYMER-METAI.  INTERFACE 
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Atoms  not  in  contact  with  the  source  (bottom  of  Fig.  7)  were  removed  from  the  electron 
micrograph  to  reveal  the  highly  ramified  nature  of  the  interface.  The  fractal  dimension  of 
the  thinly  sliced  interface  diffusion  front  was  determined  to  be  about  7/4  in  agreement  with 
Sapoval's  theory.  Details  of  this  analysis  are  given  in  references  34  and  36. 

The  ability  to  control  the  interface  roughness  and  molecular  connectivity  has 
interesting  implications  for  the  mechanical  and  electrical  property  development  of 
interfaces.  Roughness  promotes  mechanical  interlocking  but  fractal  characteristics  with 
holes  reduce  electrical  conductance  and  affect  signal  speed  in  the  metal  layer.  The 
thickness  of  the  conducting  strip  is  determined  by  the  position  of  the  fractal  diffusion  front. 
The  fractal  analysis  of  interfaces  has  been  found  useful  by  us  for  the  examination  of 
structure  development  at  reacting  polymer-cement  interfaces  and  for  investigatin 
biodegradation  of  composites  containing  degradable  and  non  degradable  components,  e.g., 
starch  and  polyethylene  blends.  Biodegradation  of  such  blends  involves  a  percolation 
invasion  process  for  materials  in  which  the  starch  molecules  are  uniformly  distributed  in 
the  PE  matrix.  However,  when  flow  induced  concentration  profiles  develop,  the  starch 
near  the  surface  is  removed  by  microbes  while  the  remaining  material  is  encapsulated 
within  the  PE  matrix.  The  surface  separating  the  accessed  material  from  the  encapsulated 
material  is  similar  to  the  fractal  diffusion  front  shown  in  Fig.  7. 


Computer  Analysis  of  Polymer-Polymer  Interdiffusion 

A  computer  simulation  of  reptating  random-coil  chains  diffusing  on  a  square  lattice 
across  an  interface  plane  was  conducted  (34).  An  IBM-AT  microcomputer  attached  to  a 
Cray  XMP-48  was  used  to  simulate  the  Brownian  motion  of  the  interpenetrating  chains. 
The  object  of  this  study  was  to  (a)  investigate  the  scaling  laws  derived  for  polymer 
interdiffu  Jon  shown  in  Table  1,  (b)  evaluate  the  non-Fickian  concentration  profiles  for 
polymer  diffusion  at  times  less  than  the  reptation  time,  and  (c)  to  examine  fractal 
characteristics  of  polymer  diffusion  fronts. 

Brownian  motion  was  accomplished  by  an  algorithm  which  randomly  chose  a  chain  end 
and  added  a  new  monomer  in  one  of  three  possible  new  directions  and  subtracted  a 
monomer  from  the  other  end.  The  process  was  then  repeated  for  each  chain.  The  chain 
end  distribution  function  was  chosen  to  be  uniform  behind  the  interface  plane.  The  basic 
reptation  dynamic  results  were  checked  for  consistency;  namely,  that  (1)  the  relaxation 
time  behaved  as  Tr  ~  M3,  (2)  the  center-of-mass  diffusion  coefficient  behaved  as  D  =  M- 
2,  and  (3)  that  the  monomer  displacement  r(t)  behaved  as, 


{ <r(t)  -  r(0)>2}n  *  tl/2. 

(t<Tr) 

(29) 

{ <r(t)  -  r(0)>2}n  *  t 

0>Tr) 

(30) 

where  r(t)  is  the  position  vector  of  the  nth  monomer  at  time  t.  The  monomer  displacement 
average  is  determined  over  all  monomers  on  each  chain. 

The  chain  density  was  chosen  to  be  inversely  proportional  to  the  molecular  weight  and 
the  monomer  density  in  the  virgin  state  was  chosen  to  consist  of  an  average  of  five 
monomers  per  lattice  site  to  account  for  interpenetration  of  random-coil  chain  statistical 
segments. 

At  t  =  0,  the  reptating  chains  were  allowed  to  diffuse  across  the  interface  plane  and 
formed  ramified  structures  as  shown  in  Fig.  7.  The  polymer  interface  is  seen  to  consist  of 
features  resembling  lakes,  islands  and  seashore  similar  to  the  previous  atomic  diffusion 
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case.  However,  many  of  the  holes  near  the  x  =  0  plane  are  not  due  to  the  fractal  nature  of 
diffusion  but  rather  are  due  to  the  monomer  correlated  motion  effect  producing  disconti¬ 
nuities  in  the  interface  as  shown  in  Fig.  5.  The  computer  simulation  gave  further  support  to 
the  scaling  laws  shown  in  Table  1  and  the  concentration  profile  proposed  by  Zhang  and 
Wool,  Eq.  21.  Fractal  properties  of  polymer  interfaces  are  discussed  in  the  next  section. 


Fractal  Polymer  DifTusion  Fronts 

The  analysis  of  the  fractal  nature  of  polymer  diffusion  fronts  (34)  is  complicated  by  the 
interpenetrated  nature  of  the  random-coil  chains  allowing  several  statistical  segments  or 
monomers  to  occupy  a  single  lattice  site.  Fig.  7  shows  a  typical  ramified  diffusion  field  for 
the  simulated  interface.  The  diffusion  front  was  obtained  using  the  same  technique 
employed  by  Sapoval,  et  al  (30).  The  front  is  potentially  more  ramified  initially  due  to  the 
connectivity  of  monomers  within  a  given  chain.  Fig.  8  shows  the  number  of  monomer  on 
the  diffusion  front,  Nm,  as  a  function  of  the  average  inter  diffusion  distance,  <X>,  for 
random-coil  chains  of  molecular  weights  ranging  from  60  to  100  statistical  units.  Two 
regions  of  behavior  are  observed. 

For  t  <  Tr,  the  front  is  discontinuous  for  reasons  discussed  above  and  Nm  can  be 
approximated  by  summing  average  contributions  from  single  chains  Vl(t)  as, 

Nm  -  n(t)Vl(t),  (31) 

where  n(t)  is  the  number  of  chains  crossing  unit  area  (or  unit  length)  of  the  initial  interface 
plane.  From  Table  1,  it  can  be  seen  that  n(t)  scales  as  n(t)  ~  X(t)/M  where  X(t)  is  the 
average  monomer  diffusion  distance  (equivalent  to  <X>  used  in  Fig.  8).  The  contribution 
to  the  diffusion  front  from  a  single  minor  chain  is  given  by  the  square  root  of  the  contour 
length  which  is  proportional  to  X(t).  Substituting  in  Eq.  31  we  obtain, 

Nm  -  X(t)7M.  (t  <  Tr)  (32) 


Fig.  8.  The  number  of  monomers  on  a  polymer  interface  diffusion  front  is  plotted  vs.  the 
average  monomer  interpenetration  depth,  <x> ,  for  molecular  lengths  ranging  from  60  to 
100  times  the  critical  entanglement  length.  The  reptation  algorithm  was  used  to  simulate 
Brownian  motion  of  linear  flexible  chains  (Wool  and  Long). 
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Fig.  8  shows  that  the  plot  of  log  Nm  vs.  log  X(t)  has  a  slope  of  2  for  all  molecular 
weights  and  that  at  constant  X(t),  Nm  varies  approximately  as  the  inverse  of  M.  However, 
for  times  less  than  the  reptation  time,  the  diffusion  front  is  not  considered  to  be  fractal 
since  the  criterion  of  self-similarity  is  not  obeyed  due  to  discontinuities  in  the  interface 
plane.  This  concept  may  need  to  be  modified  when  consideration  is  given  to  short  range 
segmental  motion  proposed  by  Zhang  and  Wool  (15).  Rouse-like  motion  of  short  chain 
segments  was  not  considered  m  this  simulation  and  the  results  reflect  the  concentration 
profile  for  monomers  on  the  primitive  path  of  the  reptating  chain.  The  time-dependent 
scaling  law  for  Nm  at  t  <  Tr  is  seen  from  Table  1  to  be  identical  to  that  for  bridges,  p(t). 

At  t  =  Tr,  since  Xqo2  «  M,  Nm(Tr)  from  Eq.  32  becomes  independent  of  molecular 
weight.  When  t  >  Tr,  the  number  of  monomers  on  the  diffusion  front  is  predicted  from 
Sapoval's  theory  to  behave  as, 

N,„  -  X0-43,  (33) 

and  is  shown  by  the  convergence  to  a  slope  of  0.43  for  all  molecular  weights  in  Fig.  8.  The 
number  of  lattice  sites,  Nl,  on  the  polymer  diffusion  front  was  found  to  behave  similarly  to 
the  number  of  monomers,  Nm.  A  mass-to-radius  analysis  of  the  diffusion  front  geometry 
gave  a  fractal  dimension  of  approximately  1.7.  These  results  are  in  agreement  with  the 
atomic  diffusion  fractal  analysis  in  2  dimensions.  Thus,  the  diffusion  front  for  polymers 
becomes  fractal  at  times  greater  than  the  reptation  time  and  at  diffusion  distances  greater 
than  the  end-to-end  vector  of  the  chains. 

MICROSTRUCTURAL  FRACTURE  CRITERIA 


Chain  Pullout 

The  Griffith  approach  to  brittle  fracture  is  given  by  (37) 

(dU/da)^  >  2r,  (34) 

where  U  is  the  stored  strain  energy  consumed  at  constant  displacement,  6,  to  advance  a 
crack  by  the  increment,  da,  and  create  surface  energy,  2rda.  The  fracture  energy,  or  the 
critical  strain  energy  release  rate,  Gjc,  is  related  to  T  by 

GIc  =  2E  (35) 

In  principle,  if  one  independently  knows  the  surface  energy,  T,  of  a  solid  one  should  be 
able  to  predict  its  fracture  energy  using  Eq.  35.  However,  for  high  molecular  weight 
entangled  polymers  the  experimental  values  are  typically  of  the  order  Gjc  *»  1000  J/m2  and 
T  ~  10"^  J/m2.  Closer  agreement  with  the  Griffith  theory  is  obtained  for  brittle  fracture  of 
very  low  molecular  weight  unentangled  polymers  (38,  39). 

The  Griffith  idea  has  considerable  value  if  one  can  determine  how  stored  energy  is 
consumed  in  forming  the  fracture  surface.  To  address  welding  and  healing  problems  a 
stored  strain  energy  approach  to  fracture  was  adopted  which  considers  both  chain  pullout 
via  disentanglement  and  chain  fracture  mechanisms  at  the  interface.  Consider  a  volume 
element  containing  a  small  area  of  interface  as  shown  in  Fig.  9  subjected  to  normal  stresses, 
ajj,  and  shear  stresses,  ry.  The  strain  energy  density,  U,  of  the  element  is  given  by, 

U  =  (o\x  +  a2yy  +  r;2zz)/2E 
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'  (°xx°yy  +  °yy°zz  +  °xx°zz)  V/E  (36) 

+  (T2xy“*"T2xz+T2yz)/2/*s> 

where  E,  v  and  /rs  are  Young's  modulus,  Poisson's  ratio  and  shear  modulus,  respectively. 

a 


Fig.  9.  A  segment  of  a  stressed  chain  at  an  interface  is  shown. 

We  first  consider  a  uniaxial  stress  field  with  a  =  a **  and  all  other  stress  terms  equal  to 
zero.  In  that  case,  the  strain  energy  density  is 

U  =  a2/2E.  (37) 

We  can  reasonably  assume  that  the  modulus  behaves  as,  (1,  40) 

E  -  t°M°;  (38) 

Ex  *  M°,  (39) 

such  that  E  can  be  assumed  to  remain  constant  during  healing  and  fracture. 

The  number  of  chains  per  unit  volume,  Nv,  is 


Nv  =  o  Na/M,  (40) 

where  o  is  density  and  Na  is  Avogadro's  number.  The  strain  energy  per  chain,  Uc,  is 
therefore  obtained  as, 

Uc  =  U/Nv  =  a2M/(2E  o  Na).  (41) 

The  strain  energy  associated  with  a  segment  of  chain  of  length  1,  is  similarly  derived  as 

Uc  -  a2l.  (42) 

The  stored  elastic  energy  can  be  used  to  either  pull  the  chain  out  from  the  interface  or 
fracture  the  chain.  A  fracture  criterion  for  chain  pullout  can  be  written  as 

Uc  >  Up  (43) 

where  Up  is  the  energy  to  pull  a  chain  out  of  its  tube  with  a  force  f,  and  velocity,  dl/dt. 
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The  pullout  force  is  given  by  f  =  /i  dl/dt,  where  /t  is  the  friction  coefficient  for  the  chain 
segment  of  length,  1.  The  friction  coefficient  for  the  segment  is  related  to  the  monomer 
friction  coefficient,  /io,  via  //  =  /rol  such  that  the  force  may  be  written  as, 


f  =  //0  1  dl/dt. 


(44) 


The  energy  required  to  pull  a  segment  out  at  constant  velocity  is  determined  from 


Up  -  Mo 


1  =  L 

l(dl/dt)  dl, 
1  =  0 


(45) 


as,  (41) 


Up  =  Mo  (dl/dt)  172,  (46) 

or 

Up  -  E 

A  similar  result  is  obtained  if  the  chain  is  pulled  out  at  constant  force  rather  than  at 
constant  velocity,  V=dl/dt.  Eq.  46  has  also  been  derived  by  Prentice  (42)  and  Evans  (43). 

Substituting  Eqs.  42  and  46  in  Eq.  43,  the  critical  fracture  stress  for  chain 
disentanglement  (and  tack)  is  obtained  as, 


a  -  (^0VEb  o  Na/Mc)V2  l(t)1/2, 
or, 


(47) 


-  X(t), 


(48) 


such  that 


a  ~  \d(t) 

where  X(t)  -  ^l(t)  is  the  average  monomer  segment  interpenetration  distance,  M0  is  the 
monomer  molecular  weight,  and  b  is  the  bond  length.  Examining  the  scaling  laws  in  Table 
1,  the  fracture  stress  or  tack  in  uniaxial  tension  should  behave  as, 


o  =  t1/4  M'1/4, 


(49) 


and  the  virgin  (green)  strength  as 
(Too  ~  M1/2. 


(50) 
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The  velocity  dependence,  V  =  dl/dt,  of  tack  and  green  strength  during  peel  adhesion 
tests  from  Eq.  47  should  behave  as, 

a  ~  V1/2.  (51) 

Fracture  by  the  chain  pullout  mechanism  is  favored  at  low  molecular  weight  (small 
l(t)),  low  deformation  rates  (small  V),  high  temperature  (low  fi0),  short  healing  times 
(small  l(t))  and  plasticizing  agents  (low  piQ). 


Fracture  Mechanics  of  Welding 

The  Dugdale  model  (44)  was  used  as  a  first  approximation  to  determine  the  molecular 
aspects  of  the  fracture  mechanics.  In  this  model,  shown  in  Fig.  10,  a  crack  of  initial  length 
ao  propagates  through  a  plastic  line  zone  or  craze  of  length  rp,  ahead  of  the  crack  tip.  The 
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Fig.  10.  The  Dugdale  model  for  a  crack  propagating  through  a  welded  interface  of 
thickness  X  -l1'2  is  shown.  The  stress,  a ,  acting  on  the  plastic  zone  of  length  rp,  draws  the 
interface  into  craze-like  fibrils  of  constant  draw  ratio,  \,  which  break  down  to  produce  the 
crack  opening  displacement,  6  (R.  P.  Wool). 

plane  stress  field  near  the  crack  for  a  Cartesian  volume  element,  xyz,  and  angle  a  at 
distance  r  is  given  by  Rice  (45)  as, 

Oxx  =  Kj/v^jit  cos  a/ 2  (1-sin  a/2  sin  3a/2)  + . ; 

Oyy  =  Kj/V2jrr  cos  a/2  (1  +  sin  a/2  sin  3a/2)  + . ;  (52) 

°xy  =  Ki/V2jrr  cos  a/2  + . , 

where  Kj  is  the  stress  intensity  factor.  The  length  of  the  plastic  zone  is  determined  by 


rp  =  n/%  Kic2/o2c,  (53) 

where  ac  is  the  craze  or  yield  stress. 

The  critical  crack  opening  displacement  is  given  by 


<5  =  Kic2/oc  E, 


where  E  is  the  tensile  modulus. 


(54) 
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Using  the  approximation  that  ac  remains  constant  within  the  plastic  zone,  the  critical 
strain  energy  release  rate,  Gjc,  is  determined  by 


Gic  =  ocd, 


(55) 


such  that  the  fracture  energy  is  determined  as  a  critical  force  acting  over  a  critical  distance. 


Fig.  11  shows  the  essential  features  of  the  Dugdale  model  for  a  crack  propagating 
through  a  line  deformation  zone  in  styrene-isoprene-styrene  (SIS)  block  copolymer  (46). 
The  microstructure  elucidated  by  Transmission  Electron  Microscopy  (TEM)  consists  of 
styrene  spheres  (white^  with  a  diameter  of  ca  250  A  in  an  isoprene  matrix  (stained  black 
with  osmium  tetroxide).  For  the  craze  zone,  rp  ~  3  fim  and  o  =  0.2  fim.  The  craze 
propagates  in  the  stress  field  applied  normal  to  the  craze  by  incorporating  a  single  styrene 
sphere  into  the  tip  of  the  zone.  As  the  zone  propagates  it  thickens  from  250  A  to  2,000  A 
by  drawing  further  material  into  the  zone  and  elongates  by  extending  the  material  already 
in  the  zone.  At  the  crack  tip,  many  styrene  spheres  have  been  drawn  into  the  zone  and 
deformed  nonuniformly  such  that  the  crack  tends  to  propagate  through  the  center  of  the 
zone.  Similar  mechanisms  have  been  identified  for  crazes  propagating  through  pure 
amorphous  polystyrene  (47,  48).  Thus,  a  critical  stress  nucleates  the  line  zone  at  distance 
rp  from  the  crack  tip  and  the  crack  opening  displacement  is  determined  by  the  breakdown 
or  the  material  drawn  into  the  zone. 


Fig.  11.  (a)  Transmission  Electron 
micrograph  of  a  plastic  zone,  rp  =  3  fim, 
in  a  mechanically  stressed  Styrene-Iso 
prene-Styrene  thin  film.  The  white 
spheres  (styrene)  have  a  diameter  of 
-250  A.  The  number  of  spheres  in  the 
zone  increases  towards  the  crack  tip 
(d=0.2  fim )  due  to  zone  thickening,  (b) 
Shows  details  at  a  crack  tip,  6=  1.0  fim, 
and  plastic  deformation  zone  extending 
into  the  undeformed  S-I-S  sample. 
(Dolman,  Robertson  and  Wool). 


To  evaluate  the  molecular  aspects  of  the  welding  problem  we  first  calculate  oc  using 
the  strain  energy  approach  such  that  the  stored  energy  per  molecule  in  an  xy-plane  stress 
field  is, 


Uc  =  (o2xx/2E  +  a2 yy/2E  +  o2x y/2/*)/Nv.  (56) 

Substituting  for  the  stress  terms  from  Eq.  52  in  Eq.  56,  Uc  is  found  to  be  dominated  by  the 
Oyy  term  and  can  be  well  approximated  by 

uc  =  °2yy  KO,  (57) 

where  l(t)  is  again  the  average  interpenetration  length.  Using  the  same  criterion  for  chain 
pullout,  Eq.  46,  and  letting  Oyy  =  ac  at  critical  conditions,  then 

ac  -  1  (t)V2,  (58) 

which  is  similar  to  the  result  obtained  for  the  uniaxial  case.  For  glassy  polymers,  ac  can  be 
related  to  the  craze  stress  which  differs  from  the  shear  yielding  stress,  Oy. 

The  critical  crack  opening  displacement,  (5,  is  determined  by  the  breakdown  of  the 
entanglement  network  in  the  deformation  zone  and  the  extent  to  which  the  zone  can 
thicken  as  depicted  in  Fig.  10.  During  welding,  the  zone  consists  of  interdiffused  short 
chain  segments  mixed  with  longer  chains  with  nonequilibrium  configurations  and  the 
entangled  structure  is  complex.  The  crack  opening  displacement  consists  of  two 
components:  a  component  from  the  virgin  material,  (3V,  and  a  component  from  the  the 
diffuse  interface,  (5j,  such  that  6  =  <5V  +  dj.  At  low  welding  strengths  the  contribution  from 
the  weaker  interface  dominates  since  the  energy  required  to  thicken  and  propagate  the 
zone  into  the  virgin  material  is  too  high.  The  crack  opening  displacement  is  then 
determined  by  the  breakdown  of  the  entanglement  network  in  the  interface  by  the 
following  approximation  (50), 

6  ~  tffl(t)1/2,  (59) 

where  ctf  is  the  constant  uniaxial  draw  ratio  to  tighten  the  slack  between  entanglements. 

At  this  point  we  briefly  consider  an  entanglement  model  which  can  be  used  to  visualize 
the  breakdown  process  of  the  deformation  zone  at  the  crack  tip  shown  in  Fig.  10.  An 
entangled  amorphous  linear  chain  network  is  shown  schematically  in  Fig.  12.  The  bridge 
theory  of  connectivity  in  an  amorphous  network  proposed  by  Wool  (3, 4)  requires  that  the 
number  of  bridges  p,  crossing  any  load  bearing  plane  exceeds  the  number  of  chains  n, 
intersecting  the  plane.  In  terms  of  a  single  chain  argument,  this  implies  that  each  chain 
must  be  long  enough  to  create  one  link  or  bridge,  i.e.,  the  chain  minimally  forms  one  circle 
or  knot  with  its  surroundings  as  shown  in  Fig.  12.  The  critical  entanglement  molecular 
weight  Mc,  occurs  at  p  =  n  and  is  given  by  (3,  4) 


Mc  =  30.89(zb/c)2Mo  CMj  (60) 

in  which  z,b,c,Mp,  Coo  and  j  are  the  number  of  monomers  per  c-axis  length  (see  examples 
below),  the  bond  length,  the  c-axis  length,  the  monomer  molecular  weight,  the 
characteristic  ratio  and  the  number  of  bonds  per  monomer,  respectively.  For  example,  with 
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polystyrene,  b  =  1.54  A,  c  =  6.5  A  for  an  isotactic  chain,  z  =  3  monomers  per  helical 
repeat  distance,  M0  =  104,  Coo  =  10,  j  =  2  c-c  bonds  per  monomer,  such  that  Eq.  60 
predicts  Me  =  32,000.  This  compares  well  with  the  critical  entanglement  molecular  weight 
for  the  cross-ovenof  the  zero  shear  melt  viscosity  from  Rouse  (rj  ~  M)  to  entangled 
behavior  (t]  ~  M3*4).  For  polyethylene,  b  =  1.54  A,  c  =  2.55  A,  z  =  1,  M0  =  28,  Cco  = 

6.7,  and  j  =  2  such  that  Mg  =  4,000.  The  calculated  values  compare  well  with  experiment 
and  can  be  approximated  for  most  vinyl  polymers  by  Mc  ~  30  Coo  M0. 

Applying  this  model  to  examine  disentanglement  as  a  function  of  draw  ratio,  a,  the 
number  of  chains  crossing  any  plane  normal  to  the  applied  strain  increases  as,n(a)  =  an  ^ 
and  the  number  of  bridges  increases  more  slowly  as  p(a)  =  pma/(a2+  \/a)l'A.  At  the 
critical  point  of  p(a)  =  n(«),  the  critical  draw  ratio  is  determined  by  ac2  +  1  /ac  =  2M/MC 
such  that, 

ac  =  (8/3  M/Mc)1/2  cos{l/3  cos'1  [-1/2(3MC/2M)3/2]}  (61) 

A  useful  approximate  solution  to  Eq.  61  at  M  »  Mc  gives  the  critical  draw  ratio  for 
disentanglement  as 


«c  =  (2M/Mc)V2  (62) 

In  a  welding  interface,  the  interdiffused  contour  length  provides  new  entanglements  which 
decrease  with  draw  ratio  such  that  a  critically  connected  state  exists  at, 


«c  -  PIM/Lc)1/2  (63) 

where  Lc  is  the  entanglement  contour  length  corresponding  to  Mc. 

Entanglement  models  (49,  50)  predict  that  the  strain  hardening  fibril  draw  ratio  is 
constant  with  values  of  =  4  which  is  in  agreement  with  TEM  studies  of  craze  fibrils  by 
Donald  and  Kramer  (49).  When  ac<af,  Rouse-like  retraction  of  the  stretched 
entanglements  will  relax  the  contour  length  such  that  the  critically  connected  state 
predicted  by  Eq.  63  exists.  An  instability  arises  which  results  in  the  completion  of  the 
fracture  process.  Substituting  for  a  and  6,  the  critical  strain  energy  release  rate  is 
determined  by 
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GIc  —  t1/2  (64) 

Goo  ~  M.  (M  <  M*)  (65) 

The  critical  stress  intensity  factor,  Kjc,  is  related  to  the  fracture  energy  via 

Gk  =  Kdc2/E,  (66) 

and  is  expected  to  increase  with  welding  time  as 

KIc  -  t1/4  M'1/4;  (67) 

Kkoo  -  M]/2.  (68) 


If  the  chains  in  the  deformation  zone  cannot  disentangle  completely  and  bond  rupture 
occurs  then  the  scaling  laws  become  more  complicated.  For  welding  of  glassy  polymers, 
fracture  involves  a  mixture  of  disentanglement  and  bond  rupture  mechanisms,  de  Gennes 
(26)  has  suggested  that  molecular  bridges  could  control  the  strength,  which  is  similar  to  the 
crossing  density  suggested  by  Prager  and  Tirreli,  (7)  and  which  is  similar  to  the  total 
number  of  entanglements  per  chain  per  unit  area  suggested  by  Kausch  (52).  In  that  case, 
the  fracture  energy  is  given  as 

GIc  =  t]/2  M'3/2;  (69) 

GjCoo  ~  (70) 

The  above  relations  for  welding  and  fracture  will  be  explored  in  the  experimental  section. 

Stages  of  Healing 

Healing  of  polymer  interfaces  has  been  described  by  Wool  and  O'Connor  in  terms  of 
several  stages  involving  (1),  (a)  surface  rearrangement,  (b)  surface  approach,  (c)  wetting, 
(d)  diffusion  and  (e)  randomization.  Most  of  the  above  discussion  focussed  on  the 
influence  of  diffusion  on  mechanical  properties.  However,  several  other  factors  can  affect 
the  welding  process  and  are  discussed  below  as  the  "stages  of  healing". 

(a)  Before  the  surfaces  are  contacted,  one  should  consider  the  roughness  or 
topography  of  the  surface  and  how  it  changes  with  time,  temperature,  and  pressure 
following  contact.  In  fractured  polymers,  rearrangement  of  fibrillar  morphology,  etc.,  can 
affect  the  rate  of  crack  healing.  Chain  end  distributions  near  the  surface  can  change  as 
molecules  diffuse  back  into  the  bulk  polymer.  Spatial  changes  of  the  molecular  weight 
distribution  may  also  occur,  e.g.,  where  the  low  molecular  weight  species  preferentially 
migrate  to  the  surface.  Chemical  reactions,  e.g.,  oxidation  and  crosslinking  can  occur  on 
the  surface  and  complicate  the  molecular  dynamics  of  diffusion.  Many  other  processes, 
e.g.,  molecular  orientation  changes,  can  contribute  to  the  stage  of  surface  rearrangement. 
Each  material  and  experimental  technique  usually  possesses  unique  surface  rearrangement 
processes  which  may  need  to  be  quantified.  O'Connor  (57)  found  significant  effects  of 
surface  rearrangement  on  the  healing  rate  of  lightly  cross-linked  polybutadiene  such  that 
the  fracture  energy  at  constant  healing  time  decreased  with  increasing  surface 
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rearrangement  time  (prior  to  contact).  This  effect  can  be  explained  by  the  chain  ends 
diffusing  back  into  the  bulk  such  that  they  were  not  as  readily  available  for  diffusion  during 
healing.  Recent  SIMS  analysis  by  McGarel  and  Wool  (unpublished)  of  low  molecular 
weight  deuterated  polystyrene  chains  on  the  surface  of  high  molecular  weight  polystyrene 
support  this  concept.  They  found  that  the  low  molecular  weight  species  will  not 
preferentially  migrate  to  the  surface  but  wili  distribute  uniformly  in  the  bulk  material. 

(b)  Surface  approach  considers  time-dependent  contact  of  the  different  parts  of  the 
surfaces  to  create  the  interface.  For  example,  in  crack,  (1, 18)  craze,  (53,  54)  and  void  (55, 
56)  healing,  contact  may  be  achieved  at  different  locations  at  different  times  in  the 
interface  depending  on  the  closure  mode  (1).  For  example,  slow  closure  of  a  double 
cantilever  beam  crack  would  result  in  different  extents  of  healing  along  the  closed  crack. 
This  stage  typically  contributes  as  a  boundary  value  problem  to  the  other  stages  of  wetting 
and  difft  sion. 

(c)  Wetting  can  occur  in  a  time-dependent  fashion  at  the  interface.  For  our  purposes 
we  provide  a  brief  phenomenological  description  of  wetting  to  illustrate  potential  problems 
in  evaluating  the  time-dependence  of  welding.  Fig.  12  shows  a  schematic  region  of  the 
plane  of  contact  of  a  polymer  interface  (1).  Due  to  surface  roughness,  etc.,  good  contac* 
and  wetting  are  not  achieved  instantaneously  at  all  locations.  Typically,  wetted  "pools"  are 
nucleated  at  random  locations  at  the  interface  and  propagate  radially  until  coalescence  and 
complete  wetting  are  obtained.  This  problem  has  been  treated  phenomenologically  as  a 
two-dimensional  nucleaiion  and  growth  process  such  that  the  fractional  wetted  area,  <I>(t), 
is  given  as  (1), 


<I>(t)  =  1  -  exp(-ktm)  (71) 

where  k  and  m  are  constants  depending  on  the  nucleation  function  and  radial  spreading 
rates. 


Contact  theories  proposed  by  Anand  (72)  and  others  suggest  that  complete  strength 
may  be  obtained  when  the  interface  has  wetted  at  O  =  1.  Others  argue  that  interdiffusion 
is  necessary  for  strength  development.  The  time  dependence  of  viscous  flow  to  promote 
contact  and  that  of  interdiffusion  may  be  comparable  since  they  are  subject  to  the  same 
molecular  dynamic  processes.  We  will  demonstrate  in  later  sections  using  both  tack 
experiments  and  incompatible  interface  weld  strength  data  that  diffusion  plays  a  very 
important  role. 

(d)  The  diffusion  stage  has  been  discussed  above  with  respect  to  the  instantaneous 
wetting  condition.  However,  in  the  presence  of  a  time  dependent  wet.ing  function,  Eq.  71, 
we  see  from  Fig.  13  that  diffusion  will  have  progressed  to  different  extents  in  different 
areas  of  the  interface.  If  the  intrinsic  diffusion  function,  H(t)  as  given  by  Eq.  23,  does  not 
change  its  nature  with  time  due  to  the  other  stages,  then  the  net  diffusion,  H'(t)  can  be 
expressed  as  the  convolution  product. 


H'(t)  = 


(t-r)  d<I>(r)/dr  dr. 


(72) 


j 


-OC 


where  r  is  a  "dummy"  variable  on  the  time  axis.  The  convolution  process  for  typical  <F(t) 
functions  may  mask  the  time-dependence  of  the  intrinsic  diffusion  function  and  related 
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mechanical  properties.  This  part  of  the  problem  can  be  solved  mathematically  by  letting 
the  wetting  function  be  a  Dirac  delta  function, 


d  <I>(t)  =  6  (r), 
dt 


(73) 


such  that  Eq.  72  reduces  to  H'(t)  -  irl(t).  This  can  be  attempted  experimentally  by 
obtaining  instantaneous  wetting  of  atomically  smooth  surfaces  under  moderate  contact 
pressure  at  temperatures  above  Tg.  However,  we  have  found  that  the  wetting  problem  can 
result  in  difficu'ty  with  the  analysis  of  most  short  time  welding  data. 


Fig.  13.  A  schematic  view  of  the  stages  of  healing  in  an  interface.  Contacted  regions  of 
radius  r  grow  and  coalesce.  Interdiffusion  occurs  only  in  the  wetted  regions  while  surface 
rearrangement  occurs  in  the  non-wetted  areas.  (Wool  and  O'Connor) 

(e)  The  randomization  stage  refers  to  the  equilibriation  of  the  nonequilibrium 
conformations  of  the  chains  near  the  surfaces  and  in  the  case  of  crack  healing  and 
processing,  the  restoration  of  the  molecular  weight  distribution  and  random  orientation  of 
chain  segments  near  the  interface.  The  conformational  relaxation  is  of  particular 
importance  in  the  strength  development  at  incompatible  interfaces. 

Tne  stages  of  crack  healing  can  have  interactive  time-dependent  functions  such  that  the 
welding  and  tack  problems  can  consist  of  processes  involving  five-way  interconvoluted 
functions  resulting  in  mechanical  properties  whose  time  dependence  may  not  be  readily 
interpretable.  Thus,  great  care  must  be  taken  in  conducting  welding  experiments  which  are 
designed  to  critically  explore  molecular  theories  of  strength  development. 

HEALING  EXPERIMENTS 


Tack  and  Green  Strength 

The  critical  experiment  is,  to  show  that  the  tack  or  stress  at  fracture  of  uncured  linear 
elastomers  behaves  as  a  ~  t*/4M '  /4V1/z!  and  that  the  green  or  virgin  strength  behaves  as 
Ocn  ~  (MV)^/  .  In  addition,  the  self-diffusion  coefficient  should  be  measurable  as  D  = 
A/M2  where  A  «s  a  constant  derivable  from  mechanical  data.  The  concurrence  of 
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predictions  derives  from  the  molecular  dynamics  analysis  of  tack  and  green  strength.  These 
predictions  apply  to  all  amorphous  linear  polymers  whose  chain  configurations  can  be 
described  by  Gaussian  statistics  and  whose  fracture  behavior  is  dominated  by  chain 
disentanglement  mechanisms. 

Fig.  14  shows  tack  data  obtained  by  Skewis  (61)  for  several  polymer-polymer  pairs 
where  we  have  replotted  the  tack  (units  of  force  applied  to  a  cross-sectional  area  of  0.4 
cm2)  vs.  t1/  .  The  original  data  were  nonlinear  when  plotted  on  a  t1/^  scale  to  investigate 
diffusion  mechanisms.  However,  the  fracture  energy  behaves  as  Gv  ~  t ^  2  while  the  tack 
behaves  as  a  «  t1/4  as  evidenced  by  the  linear  plot  in  Fig.  14.  The  §BR/Butyl  pair  are 
incompatible  and  the  long  time  Muilibrium  strength  is  considerably  weaker  than  either  of 
the  symmetric  pairs.  Similar  c1/4  results  for  tack  have  been  reported  by  Wool  and 
O'Connor  (1,  40)  for  polybutadiene,  and  by  Voyutskii  (11)  for  polyisobutylene. 

CONTACT  TIME  (SEC) 


Fig.  14.  The  tack  or  breaking  load  vs.  t^/4  for  several  polymer-polymer  uncured  pairs  as 
investigated  by  Skewis  (61).  Tack  measurements  were  done  with  a  contact  load  of  1000  g 
and  a  contact  area  of  0.40  cm.  The  viscosity-average  molecular  weight,  Mv,  of  the  SBR 
was  260,000  and  225,000  for  the  butyl  rubber. 

The  molecular  weight  dependence  of  tack  and  green  strength  was  investigated  by 
Forbes  and  McLeod  (62)  for  fractionated  samples  of  Hevea.  Their  results  are  shown  in 
Fig.  15  in  which  the  upperline  represents  the  green  strength  and  the  bottom  line  represents 
the  tack  stress  evaluated  at  a  constant  contact  time  of  t,-  =  30  sec  for  each  molecular 
weight.  The  green  strength  curve  was  obtained  theoretically  using  (1) 


a oo  =  k^M1/2  -  Mjl/2),  (74) 

where  kj  =  993  Pa/M^/2  and  M;  V2  =  284.1  is  the  intercept  on  the  Ooo  =  0  axis  of  the 
plot  of  cr oo  vs.  M1/  .  The  curve  through  the  tack  data  in  Fig.  14  was  determined  by  theory 

as 

a  =  aa  +  aM'1/4,  (75) 

where  a  =  5.29  MPa  M4/4  and  o0  =  80.5x10^  Pa  is  .a  constant  dependent  on  the  initial 
wetting  conditions  of  the  interface.  When  a  =  k2t1/4  as  in  Fig.  14,  the  constant,  k2,  is 
related  to  the  constant,  a,  in  Eq.  75  by 
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k2  =  a/(M1/4  tcl/4x 


(76) 


where  in  this  case  tc  =  30  sec.  Considering  the  difficulties  in  performing  such  experiments, 
the  agreement  between  theory  and  experiment  in  Fig.  15  is  satisfactory.  If  the  strength  is 
controlled  by  bridges,  then  a  ~  M'2*/4  at  constant  tc  and  a<x>  ~  M°.  These  predictions  are 
inconsistent  with  the  data  in  Fig.  15.  It  has  been  suggested  that  the  tack  stress  is 
independent  of  molecular  weight  but  in  fact  it  behaves  approximately  as  M*1'4  jn  close 
agreement  with  disentanglement  theories.  This  is  a  weak  molecular  weight  dependence 
but  related  mechanical  properties  such  as  fracture  war k,  Gi^,  should  behave  with  a 
stronger  molecular  weight  dependence,  as  Gic  ~  t1' . 


In  Fig.  15  the  green  strength  is  seen  to  rise  steadily  up  to  M  *  1.5  x  10^  at  which  point 
it  remains  fairly  constant  and  bond  rupture  is  expected  to  dominate  the  fracture  mechanism 
by  regulating  the  extent  of  disentanglement.  Given  sufficient  time,  the  tack  stress  will 
increase  to  the  green  strength  value  by  continued  diffusion.  Using  tracer  diffusion 
methods  Skewis  (63)  measured  the  self-diffusion  coefficient  of  natural  rubber  at  25°C  as  D 
=  3x10'44  cm2/sec  using  samples  with  a  viscosity-average  molecular  weight  of  Mv  = 
234,000.  The  polydispersity  of  his  samples  was  not  known.  The  self-diffusion  coefficient 
from  healing  data  is  given  by  (1-3) 


D  =  (k2/k1)4b2Coo/6MM0,  (77) 

where  Coo  is  the  characteristic  ratio  for  the  dimensions  of  a  random-coil  chain.  Letting  M 
=  234,000,  Mp  =  68  for  polyisoprene,  b  =,4.05  A  from  the  c-axis  of  the  cis-l,4-polyiso- 
prene  unit  cell,  and  k2  =  0.103  MPa/sec1/4  from  Eq.  76.  we  obtain  D  ~  10" 14  cm2/sec 
which  is  of  the  same  order  of  magnitude  as  D  =  3  x  10’14  cm2/sec  obtained  by  Skewis. 


Fig.  15.  Tack  (circles)  and  green  strength  (triangles)  as  a  function  of  viscosity  average 
molecular  weight  for  fractionated  samples  of  natural  rubber  (Forbes  and  McLeod)  (62). 
The  solid  line  for  green  strength  was  theoretically  obtained  using  Eq.  72  and  the  dashed 
line  for  tack  using  Eq.  73.  The  green  strength  was  evaluated  at  a  uniaxial  test  speed  of 
26.67  cm/min  at  25°C,  and  the  tack  was  evaluated  at  a  constant  contact  time  of  30  sec  for 
each  sample.  Filled  and  unfilled  data  points  refer  to  unmillled  and  milled  natural  rubber, 
respectively. 

The  reptation  time,  Tr,  is  obtained  from  the  self-diffusion  coefficient  as 


Tr  =  R2/(3;r2D),  (78) 

where  R2  is  the  end-to-end  vector  given  by  R2  =  Mj/M0  Coo  b2.  For  natural  rubber  the 
healing  time  is  approximated  by 
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Tr  =  2.48  x  10-!5  sec.  (79) 

When  M  =  234,000,  Tr  =  32  sec  and  in  Fig.  15,  this  material  should  be  in  the  green  state  at 
t^  =  30  sec  as  observed.  The  healing  time  required  for  the  samples  with  M  =  1,000,000  is 
about  40  min.  These  results  differ  with  those  of  Roland  and  Bohm  (69)  who  found  that  the 
time  required  to  achieve  maximum  autohesive  strength  of  polybutadiene  was  orders  of 
magnitude  longer  than  the  time  required  for  the  chains  to  diffuse  a  distance  approximately 
equal  to  the  end-to-end  vector,  R.  Chain  branching  and  nonuniform  wetting  at  the 
interface  may  have  contributed  to  this  difference. 

Stacer  and  Schreuder-Stacer  (95)  investigated  the  time  dependence  of  aiUohesion  in 
polyisobutylene.  Molecular  weights  of  their  samples  ranged  from  1.5-21  x  l(r  with  poly- 
dispersities  of  1.1  to  3.6.  T-peel  test  configurations  were  used  to  evaluate  the  peel  energy, 
G,  as  a  function  of  time,  molecular  weight,  temperature  and  pressure.  Fig.  16  shows  a  plot 
of  log  G  vs  log  t / ax  for  PIB  with  M„,  =  1.5  x  1(P.  The  contact  time  was  reduced  with 
respect  to  the  time-temperature  shift  factor,  ax,  to  create  a  master  plot  from  data  obtained 
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Fig.  16.  Comparison  of  the  predictions  of  several  theoretical  models  with  autohesion  data 
for  PIB  (Stacer  and  Schreuder-Stacer). 

at  different  temperatures.  Several  theories  were  examined.  The  scaling  law  for  peel 
energy  is  G  ~  (t/M)1/ z  as  given  by  Wool,  et  al.,  in  Eq.  64,  predicts  a  slope  of  1/2  as  shown 
in  Fig.  16.  They  conclude  that  the  1/2  power  law  slope  is  applicable  to  describe  limited 
regions  of  the  data  but  does  not  provide  an  overall  description  of  the  autohesion  curve. 

The  time  to  achieve  equilibrium  strength  (healing  time)  too,  was  predicted  by  Wool  (3)  as, 

too  =  Coo  b^M^  /  (6AMq)  (80) 

where  A  is  the  constant  in  the  self  diffusion  coefficient  relation,  D  =  A/M2.  This  relation 
was  found  to  be  in  excellent  agreement  with  the  PIB  autohesion  data.  The  molecular 
weight  dependence  of  healing  was  not  properly  analyzed  in  this  work  since  the  M'l/Z 
dependence  of  the  healing  rate  was  incorrectly  compared  to  the  green  strength  energy,  Goo- 

The  monodisperse  and  polydisperse  chain  interdiffusion  theories  of  Prager  and  Tirrell 
are  compared  with  experiment  in  Fig.  16.  This  model  determines  the  number  of  effective 
crossings  per  unit  area  where  an  effective  crossing  involves  a  minimum  interpenetration 
distance  related  to  the  entanglement  molecular  weight  spacing  which  allows  a  chain 
segment  to  participate  in  energy  dissipative  processes.  The  agreement  is  good  for  the 
monodisperse  case  and  the  polydisperse  model  reasonably  described  similar  studies.  Both 
models  predict  a  slower  rate  of  healing  (about  1/2  time  decade)  compared  to  experiment. 
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Fig.  16  also  compares  a  wetting  theory  of  Korenevskaya,  et  al.  (961,  expressed  as  a  first 
order  kinetic  process  in  terms  of  the  fractional  peel  energy  G  =  G(t)/Gx  as 


G  +  In  (1-G)  =  -  P  t//r 

where  P  is  the  applied  pressure  and  fi  is  the  zero-shear  viscosity.  This  relation  provided 
excellent  agreement  with  the  data.  Deviations  to  significantly  lower  Gjc  values  were 
observed  at  short  times  with  the  high  molecular  weight  PIB  samples. 

Since  the  diffusion  controlled  healing  times  are  comparable  to  characteristic  relaxation 
times  of  viscosity,  this  work  cannot  readily  distinguish  between  weiting  and  diffusion 
theories.  Studies  by  Wool  and  O'Connor  (21,  57)  of  tack  and  healing  in  hydroxy  terminated 
polybutadiene  (HTPB)  indicated  that  the  mechanical  property  recovery  due  to  wetting 
occurred  in  a  few  minutes  and  was  a  small  fraction  of  Goo,  while  interdiffusion  of  the 
partially  cross-linked  chains  took  several  weeks  during  which  the  majority  of  the  fracture 
energy  was  recovered.  The  separation  of  the  relaxation  times  for  wetting  and  diffusion  in 
that  case  support  the  diffusion  model. 

A  more  convincing  case  for  the  diffusion  model  is  obtained  from  studies  of 
incompatible  amorphous  polymer  interfaces  such  as  PS/PMMA.  At  long  contact  times, 
complete  wetting  or  contact  has  been  achieved  but  the  incompatibility  only  allows  limited 
interdiffusion  with  with  correspondingly  low  fracture  energy.  The  observed  fracture 
energies  of  incompatible  interfaces  by  the  Urbana  group  are  consistent  with  the  diffusion 
model  such  that  the  fracture  energy  increases  as  the  compatibility  of  the  polymer  pairs  and 
interdiffusion  increases  at  constant  contact  area.  A  possible  flaw  in  this  argument  is  that 
the  incompatibility  influences  chemical  interactions,  e.g.,  acid-base,  at  the  contact  plane 
which  would  affect  the  wetting  theory.  However,  the  contact  plane  no  longer  exists  as  such 
due  to  the  limited  diffusion  which  has  occurred. 

Pressure  affects  diffusion  by  decreasing  the  free  volume  which  slows  the  hopping 
process  necessary  to  achieve  Brownian  motion  and  subsequently  decreases  the  self¬ 
diffusion  coefficient.  This  effect  is  more  pronounced  near  the  glass  transition  temperature 
and  is  important  in  polymer  melt  processing  where  large  hydrostatic  pressures  can  be 
encountered.  It  is  not  very  important  in  normal  tack  experiments  where  the  contact 
pressures  are  usually  much  less  than  a  kilobar.  However,  pressure  affects  the  wetting  stage 
by  promoting  better  contact.  Thus,  tack  measurements  should  be  dependent  on  pressure  to 
some  degree  while  wetting  is  occurring  but  should  be  largely  independent  of  pressure 
when  interdiffusion  is  occurring.  Expressions  for  the  pressure  dependence  or  Tr  and  D 
have  been  given  (1,  3). 

The  effect  of  strain  rate  and  peeling  velocity  on  tack  has  been  examined  by  many 
investigators.  Fig.  17  shows  data  obtained  by  Hamed  (64)  for  the  cohesive  strength  of  an 
SBR  elastomer  (FRS-146)as  a  function  of  peel  rate  using  a  T-peel  test  configuration.  The 
prediction  that  a  ~  V1/2  is  observed  to  be  in  agreement  with  the  high  temperature  data 
(where  chain  pullout  is  expected  to  dominate)  as  judged  by  the  slope  of  1/2  for  the  plot  of 
log  a  vs.  log  rate.  As  the  temperature  approaches  Tg,  -50°C,  chain  fracture  begins  to 
contribute  to  the  deformation  mechanism  and  the  \a/2  behavior  is  lost  as  predicted.  The 
prediction  that  o  ~  V1/2  is  also  supported  by  data  of  Tsuji,  et  al.  (65)  for  T-peel  tests  of 
polyisobutylene  and  by  data  of  Bhowmick,  et  al.  (66),  for  adhesive  tack  and  green  strength 
of  EPDM  rubber.  The  dependence  of  the  fracture  energy,  G,  on  the  rate  of  peel  is  also  G 
~  V1/2  since  the  extent  of  disentanglement  and  not  the  rate  determines  the  crack  opening 
displacement  for  the  chain  pullout  mechanism. 


Fig.  17.  Cohesive  tear  strength  of  SBR  as  a  function  of  peeling  rate  at  several  constant 
temperatures  (Hamed  and  Shieh)  (64).  The  dashed  line  represents  the  theoretical  slope  of 
1/2  and  was  added  for  comparison  purposes.  The  SBR  was  a  cold  emulsion  styrene- 
butadiene  random  copolymer  containing  23%  bound  styrene  and  had  a  Mooney  viscosity 
(ML-4,  100°C)  of  40. 

For  peel  adhesion,  the  pullout  mechanism  predicts  the  fracture  energy  to  behave  as 


G^CtV/M)1/2;  (81) 

Goo  ~  MV1/2.  (82) 

Several  other  theories  of  tack  and  adhesion  have  been  presented.  Bister,  et  al.  (67), 
predict  that  the  peel  energy  behaves  as 


G  *  M2/3  (Dt)1/2.  (83) 

Since  D  =  M'2  and  Tr  =  M3,  then  Eq.  82  predicts  that  G  ~  t1/2  M‘~/3  and  Goo  ~  M"1/^. 
However,  this  molecular  weight  dependence  cannot  be  supported  by  data  reported  here  or 
elsewhere.  Most  theories  predict  the  t1/2  dependence  but  as  seen  from  Table  1,  this  is  not 
a  unique  function. 

Wu,  et  al.  (68)  predict  that  the  adhesive  strength  is  determined  by  G  *  n(t)X(t)  which 
behaves  in  an  identical  manner  to  bridges  (Table  1)  as  G  *  t1/2  M'-V2  and  Goo  ~  M°. 
The  molecular  weight  dependence  of  these  relations  is  not  supported  by  experiments. 
However,  Wu,  et  al.,  also  suggest  that  the  number  of  chains,  n(t),  saturates  quickly  and  that 
G  depends  on  center-of-mass  motion  over  larger  distances  as  G  ~  t1/2  M  .  We  disagree 
with  the  latter  suggestion  since  n(t)  increases  simultaneously  with  the  strength  for  times  t  < 
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Vasenins  (70)  kinetic  theory  of  adhesion  suggests  that  the  tack  stress  behaves  as 
a  »  t1/4  M'5/3  V. 

While  the  time-dependence  is  supported  by  experiment  the  molecular  weight  and  rate 
dependence  are  at  variance  with  the  body  of  data.  Many  other  theories  of  tack  and 
adhesion  are  discussed  in  references  10  and  71.  Kausch  aiso  presents  a  review  of  crack 
healing  studies  in  reference  18. 


Summary  of  Tack  and  Green  Strength 

Strength  development  at  a  polymer-polymer  interface  was  analyzed  in  terms  of  the 
dynamic  and  static  properties  of  random-coil  chains.  Interdiffusion  of  chain  segments 
across  the  interface  was  considered  to  be  the  controlling  factor  for  tack  and  green  strength 
of  uncured  linear  elastomers.  This  concept  is  similar  to  that  proposed  by  Voyutskii  but 
differs  markedly  from  contact  theories  as  proposed  by  Anand  (72)  and  others.  In  our 
approach,  time-dependent  wetting  first  occurs  followed  by  interdiffusion.  Increasing 
contact  pressure  and  temperature  should  promote  the  establishment  of  molecular  contact 
(wetting)  at  the  interface  up  to  a  saturation  point  of  complete  wetting.  However, 
interdiffusion  is  retarded  by  increased  hydrostatic  pressure  and  enhanced  by  temperature 
in  the  usual  thermally  activated  manner.  The  effect  of  pressure  on  diffusion  is  to  reduce 
the  volume  available  (depending  on  the  compressibility)  for  segmental  motion  and 
subsequently  decrease  the  diffusion  coefficient. 

Increasing  the  tack  test  temperature  increases  the  average  interdiffusion  chain  segment 
length,  as  1(T)  *  exp  -Q/2kT,  but  decreases  the  stress  required  to  pull  the  segment  out,  as 
o  ~  expQ/4kT,  where  Q  is  the  activation  energy.  Therefore,  the  tack  evaluated  at  constant 
contact  time  will  decrease  with  increasing  temperature  for  the  interdiffnsion  controlled 
process. 

The  effect  of  molecular  weight  on  tack  at  constant  contact  time,  t^  is  to  increase  the 
tack  according  to  j  for  those  samples  whose  relaxation  time  Tr  <  tc,  and  decrease 
the  tack  as  a  «  M'1/4  for  those  molecular  weights  for  which  Tr  >  t^.  The  tack  should 
reach  a  maximum  at  a  molecular  weight  corresponding  to  Tr  =  tc.  The  latter  value  can  be 
used  to  determine  the  self-diffusion  coefficient  of  the  chains  (1).  At  small  contact  times  the 
results  might  be  complicated  by  wetting  processes.  The  position  ofLthe  maximum  in  a  plot 
of  a  vs.  M  is  relatively  insensitive  to  the  contact  time  since  Tr  *  M3  and  M  at  the 
maximum  will  consequently  increase  as  t^'  .  These  predictions  appear  to  be  in 
agreement  with  much  experimental  data  reviewed  by  Hamed  (58)  and  Rhee  (59)  but  differ 
in  some  respects  from  their  own  interpretations  of  the  same  data. 

The  above  results  are  not  unique  to  elastomers  and  many  of  the  predictions  will  be 
seen  to  describe  welding  and  fracture  of  glassy  polymers  in  the  next  section. 

WELDING  OF  POLYMER  INTERFACES 
Symmetric  Amorphous  Interfaces 

Fig.  18  shows  the  welding  method  used  to  evaluate  healing  relations  for  polymer 
interfaces  (73,  74).  Monodisperse  molecular  weight  films  (or  other  materials)  are  bonded 
to  substrates,  usually  composed  of  the  same  polydisperse  molecular  weight  material,  and 
the  film  surfaces  are  then  welded  together.  The  film  surfaces  to  be  contacted  are  first 
molded  against  highly  polished  metal  plates  and  annealed  in  vacuum  to  assure  maximum 
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smoothness.  To  approximate  the  instantaneous  wetting  condition  necessary  to  evaluate  the 
effect  of  diffusion  on  welding,  the  surfaces  are  wetted  at  a  higher  temperature  for  a  short 
time,  t «  Tp  and  subsequently  healed  for  long  times  at  a  constant  lower  temperature. 

Under  these  conditions,  the  long  time  healing  data  extrapolates  to  zero  time. 
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Fig.  18.  Wedge  cleavage  method  of  testing  strength  vs.  welding  time  in  polymers. 
Monodisperse  molecular  weight  films  are  attached  to  polymer  substrates  prior  to  thermal 
welding  (Wool  and  O'Connor). 

The  stress  distribution  and  the  fracture  mechanics  of  the  wedge  cleavage  method  have 
been  examined  (75,  76).  The  wedge  is  driven  into  the  interface  spreading  the  arms  apart 
until  the  crack  of  length  a  begins  to  grow  at  a  critical  value  of  d  as  shown  in  Fig.  18.  By 
modelling  the  specimen  as  an  elastic  double  cantilever  beam  with  each  beam  partly 
supported  by  an  elastic  foundation,  the  critical  stress  intensity  factor,  Kjc,  is  obtained  at 
constant  displacement  as  (75,  76) 


KIc  =  3/4  Eh3/2(d-2h)  f(h/a)/a2.Q,  (84) 

where  h  is  the  half-width  of  the  specimen  (see  Fig.  18)  and  d  is  the  separation  of  the 
specimen  at  the  upper  end.  The  function  f(h/a)  is  given  by 


f(h/a)  .  <U°-64h/a) _ r 

1  +  1.92h/a  +  1.22(h/a)2  +  0.39(h/a)3 


(85) 


and  the  factor  Q  is  a  correction  for  the  extra  compliance  due  to  the  presence  of  the  end- 
slot  (76),  Q  =  1  +  [(h/h')3  - 1]  (ls/a)  ,  where  Is  is  the  length  and  h'  is  the  half-width  of  the 
slot,  respectively.  The  fracture  energy  is  obtained  from  Eq.  84  via  the  relation  Gjc  = 
KIc2/E.  The  advantage  of  the  wedge  cleavage  technique  is  that  the  crack  propagates  in  a 
stable  manner  such  that  multiple  data  points  can  be  obtained  from  a  single  sample  thus 
providing  good  statistics.  It  is  possible  to  obtain  the  weld  strength  as  a  function  of  welding 


34 


time  by  alternating  crack  propagation  with  healing  treatments  such  that  each  measurement 
corresponds  to  crack  advance  into  a  previously  unfractured  portion  of  the  interface. 

Fig.  19  shows  results  obtained  by  O'Connor  and  Wool  (corrected  for  initial  wetting)  for 
welding  monodisperse  molecular  weight  polystyrene  interfaces  at  120°C.  The  best  fit  to  the 
data  gives, 


KIc  *  t1/4  M‘x,  (86) 

where  x  is  in  the  range  0.31  to  0.37.  The  latter  value  can  be  compared  to  the  chain  pullout 
prediction  of  x  =  0.25  and  the  bridge  fracture  prediction  of  x  =  0.75. 
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Fig.  19.  Kjc  vs.  t^/4  for  welding  of  polystyrene  of  varying  narrow  fraction  molecular 
weights  using  the  wedge  cleave  technique  (Wool  and  O'Connor). 

Petrovska  and  Kausch  (20)  studied  crack  healing  of  PMMA  using  compact  tension 
specimens  with  molecular  weights  of  Mw  =  88,500, 162,000, 404,000  and  727,000.  The 
samples  were  fractured  at  room  temperature  and  the  recontacted  surfaces  were  healed  at 
117°C.  They  found  that  their  data  could  be  described  by  Eq.  86  with  x  =  0.14.  The  effect 
of  bond  rupture  in  the  healing  experiments  would  be  to  lower  the  exponent  for  the 
molecular  weight  dependence,  compared  to  welding  of  unfractured  surfaces.  In  crack 
healing,  the  healing  rate  is  determined  by  the  effective  molecular  weight  Mf,  in  the 
interface  of  thickness  X.  The  number  of  broken  bonds  has  been  determined  by  Wool  and 
Rockhill  (17)  as, 


Nf  =  o  Nx  X/2  (1/Mf  -  1/M)  (87) 

where  Nf  is  the  number  of  broken  bonds  per  unit  fracture  surface  area,  Na  is  Avogadro's 
number  and  M  is  the  virgin  state  molecular  weight.  Appropriate  values  are  Nf  =  5  x 
lO^/cm2  and  X  ~  1000  A  (or  appropriate  molecular  weignt  dependence).  Solving  for  Mf 
and  substituting  Mf  for  M  in  Eq.  86,  it  is  easy  to  show  that  a  lower  exponent  compared  to 
the  virgin  state  welding  case  would  be  obtained  for  the  molecular  weight  dependence  of  the 
healing  rate. 

PMMA  crack  healing  experiments  otOiConnor  (57)  indicated  that  the  time 
dependence  of  Kjc  was  non  linear  on  a  t1/4  scale  which  may  be  due  to  the  molecular 
weight  changes  on  the  crack  surfaces  and  the  presence  of  craze  fibrillar  material.  A  higher 


slope  was  observed  at  short  healing  times  which  changed  to  a  lower  slope  at  longer  times. 
Higher  molecular  weight  samples  gave  slower  healing  rates  as  observed  by  Kausch,  et  al. 
(20)  but  the  data  scatter  in  time  and  molecular  weight  did  not  permit  a  simple  scaling  law 
analysis  using  Eq.  86. 

Fig.  20  shows  Gjc  welding  data  at  115°C  by  McGarel  and  Wool  (77,  78)  for  wedge 
cleavage  welding  of  monodisperse  molecular  weight  samples  (152,000  and  400,000)  of  poly¬ 
styrene.  The  data  is  plotted  as  Gjc  vs.  tV^M-1' ^(contour  length).  The  data  exhibits 
reasonable  agreement  with  the  contour  length  scaling  law.  Very  little  superposition  was  . 
obtained  when  the  data  was  plotted  with  respect  to  the  bridge  scaling  law,  2 
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Fig.  20.  Gjc  vs.  t^M'1/2  for  welding  of  polystyrene  of  molecular  weights  152,000  and 
400,000  (Pressure  Chemical)  at  115°C  (McGarel  and  Wool). 

In  accordance  with  the  scaling  laws  for  healing,  if  G  a  l(t)  then  the  virgin  state  should 
support  GtCoo  ~  M  for  some  range  of  molecular  weight.  Fig.  21  shows  the  fracture  energy, 
Gjc,  of  polystyrene  as  a  function  of  molecular  weight  (monodisperse  fractions)  (50,  73,  77). 
An  abrupt  change  in  fracture  behavior  is  observed  in  the  vicinity  of  Mc  ~  32,000  which 
parallels  changes  in  the  zero-shear  viscosity  behavior.  For  M  <  Mc,  the  data  obtained  by 
Robertson  (79)  using  wedge  cleavage  methods  indicate  that  the  material  is  very  fragile  with 
Gir  values  of  less  than  1  J/m2.  The  fracture  energy  in  this  region  is  well  described  by  the 
linear  equation, 


GIc  =  0.23  +  3.44x10.5M  (J/m2). 


(88) 


As  M  decreases,  the  fracture  energy  approaches  the  Griffith  limit,  Gic  =  2T  where  T  is 
the  surface  energy.  For  PS  at  room  temperature  2r  =  0.08  J/m2  which  is  smaller  than  the 
extrapolated  limit  of  0.23  J/m2  for  the  fracture  energy  in  Fig.  21.  However,  this  difference 
may  be  reasonable  when  consideration  is  given  to  the  surface  obtained  by  fracture 
compared  with  the  flat  surface  normally  used  in  surface  tensiorr measurements.  The  same 
data  can  also  be  described  by  a  power  law,  Gtc  =  4.67  x  10"2MU252  J/m2.  It  has  heep 
suggested  by  Kramer  (80)  from  an.analysis  of  the  Robertson's  data  that  Gjc  *  M1' 2  based 
on  the  random-coil  radius  <S2>  V-2  ~  M1/2  controlling  the  crack  opening  displacement  at 
constant  craze  stress.  The  result,  Gjc  «  M,  which  is  consistent  with  chain  pullout 
arguments,  provides  a  slightly  better  fit  of  the  data  in  the  molecular  weight  range  M  <  Mc. 
The  data  are  too  limited  and  the  function  differences  too  small  in  the  available  data  range 
to  confidently  determine  the  better  approach. 
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In  the  molecular  weight  range  Mc  <  M  <  M*,  where  M*  ~  8MC  «  250,000,  Gjc 
increases  from  1  to  1000  J/m2  as  shown  in  Fig.  21  and  can  be  described  approximately  by 

GIc  «  (M1/2-Mc1/2)2  (89) 

The  fracture  data  are  normalized  (by  subtraction)  to  Mc.  The  physical  meaning  of  this 
subtraction  is  that  at  Mc,  the  melt  is  critically  entangled  such  that  GIc  ~  0  and  an 
instability  develops  which  does  not  require  the  removal  of  the  average  remaining  bridge  of 
length  Lc.  When  M  »  Mc  but  less  than  M*,  then  the  relation  Gr^  =  M  is  realized.  The 
uniaxialfracture  stress,  a,  and  the  critical  stress  intensity  factor,  Kjc,  both  behave  as  a,Kjc 
=  (M1/  -Mg1/2)  in  this  range  (1).  Entanglement  connectivity  arguments  based  on  percol¬ 
ation  (84, 85)  have  been  presented  to  explain  the  role  of  Mc  in  fracture  of  glassy  polymers 
/4,  50).  When  M  >  Mc,  chain  pullout  and  disentanglement  are  dominant  mechanisms  of 
fracture  in  glassy  polymers.  The  fracture  surface  morphology  consists  of  fractured 
remnants  of  craze  fibrillar  material  and  considerable  energy  is  expended  in  the  formation 
and  disentanglement  of  the  craze  microfibrils.  The  critical  energy  to  propagate  a  crack 
through  a  craze  zone  is  given  by  Gjc  =  o^d  where  oc  i>  the  stress  necessary  to  form  the 
craze  and  6  is  the  critical  crack  opening  displacement.  Both  ac  and  d  should  show,a 
dependence  on  M  and  Mc  which  behaves  approximately  as  ac  ~  d  ~  (M1/  -Mj.1/2).  Some 
support  for  this  prediction  is  provided  by  fracture  data  on  polycarbonate  obtained  by 
Pitman  and  Ward  (81). 


Molecular  Weight 

Fig.  21.  Gir  vs.  molecular  weight  for  polystyrene  in  the  virgin  state.  Data  of  Robertson 
(M  <  40,000,  filled  points)  and  O'Connor  and  Wool  (unfilled  points). 

When  M  >  M*,  the  fracture  parameters  become  independent  of  molecular  weight  as 
shown  in  Fig.  21.  Tbe  fracture  energy  attains  an  upper  bound  value  of  Gjc  =  1000  J/m2. 
Upper  bound  values  as  low  as  400  J/m2  have  also  been  reported  and  appear  to  depend  on 
the  notch  preparation  method.  If  the  notch  is  exceedingly  sharp,  the  lower  value  is 
obtained  such  that  a  single  craze  propagates  ahead  of  the  crack.  Otherwise,  craze  bundles 
and  other  dissipative  processes  can  contribute  to  the  fracture  energy.  Welding  of  polymer 
interfaces  provides  an  exceptionally  sharp  notch  by  virtue  of  the  weak  interface  prepared 
with  smooth  surfaces.  The  asymptotic  weld  strength  at  t>Tr  converges  on  the  true  upper 
bound  virgin  state  strength.  Recent  welding  results  of  McGarel  (78)  suggest  that  the  value 
of  Gjc  *=  400  J/m2  is  the  correct  upper  bound  "sharp"  notch  value  for  the  virgin  state 
fracture  energy  of  PS.  This  observation  leaves  room  for  further  interpretation  of  the 
molecular  weight  dependence  of  the  fracture  energy  depicted  in  Fig.  21. 
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The  number  of  broken  bonds  during  fracture  of  PS  in  the  uppen bound  molecular 
weight  range  was  determined  at  room  temperature  as  Nf  =  7  x  101 ' /m2  (17,  82,  83).  Bond 
rupture  contributes  about  0.0001  kJ/m2  to  the  fracture  energy  of  1  kJ/m2  and  is  essentially 
negligible.  However,  the  number  of  broken  bonds  roughly  corresponds  to  the  number  of 
bridges  crossing  the  fracture  plane  and  suggests  that  bond  rupture  controls  the  extent  of 
chain  disentanglement  and  plastic  deformation  of  the  craze  at  M  >  M*.  Thus,  in  fracture 
of  glassy  polymers,  viscous  and  plastic  deformation  processes  dominate  the  energy 
consumption  while  bond  fracture  controls  their  extent.  As  the  fracture  temperature  of  PS 
is  increased  to  T„  the  number  of  broken  bonds  decreases  (82,  83)  and  chain  pullout  is 
favored.  Complete  disentanglement  (Nf  =  0)  is  predicted  to  occur  for  PS  at  Tg  +  90°C. 

Prentice  (42)  and  Evans  (43)  have  argued  that  the  fracture  energy  should  depend  cm 
molecular  weight  via  relations  similar  to  Eq.  46  for  chain  pullout  such  lhal  Gfc  «  nooL  . 
While  they  incorrectly  assume  that  noo  ~  M°  (from  Table  1,  no,  «  M_i/Z)  they  suggest 
that  Gfc  «  Lr  «  Mz  and  show  from  PMMA  fracture  data  similar  to  that  shown  in  Fig.  21 
that  a  slope  of  2  is  obtained  in  the  range  Mc  <  M  <  M*.  A  power  law  fit  of  the  data  in  Fig. 
21  indicates  that  Gtc  «  M-3-4  but  the  slope  of  3.4  is  due  to  the  inhomogeneous  nature  of 
qic  vs.  M  and  the  data  is  more  justifiably  represented  by  relations  similar  to  Eq.  89  above. 
If  the  relation  Gic  ~  M2  were  correct,  then  one  would  expect  the  time  dependence  of 
welding  as  Gfc(t)  ~  t/M,  but  this  is  not  consistent  with  welding  data. 


Effects  of  Pressure  on  Welding 

The  effect  of  the  applied  pressure  on  welding  of  monodisperse  samples  of  polystyrene 
was  studied  by  McGarel  and  Wool  (97).  Blocks  of  52  x  52  x  3.5  mm  of  PS  were  welded  at 
different  pressures  for  30  min  at  115°C.  The  blocks  were  then  cut  into  compact  tension 
fracture  mechanics  specimens,  the  notch  was  fatigue  sharpened  at  10  Hz  to  produce  a 
starter  crack  and  the  specimens  were  fractured  at  a  rate  of  0.25  mm/min  at  room 
temperature.  When  no  pressure  was  applied  during  welding,  the  blocks  pulled  apart  due  to 
thermal  contraction  upon  cooling.  The  effect  of  pressure  on  Gic  at  constant  welding  time 
is  shown  in  Fig.  22.  The  fracture  energy  of  the  interface  is  seen  to  be  independent  of 
pressure  above  0.4  MPa  and  Gjc  values  were  in  the  range  of  327  to  385  J/m2  which  is 
below  the  virgin  fracture  value  of  700  J/m2. 

The  weld  pressure  influences  the  strength  development  in  two  ways.  First,  a  certain 
amount  of  pressure  is  necessary  to  promote  intimate  contact  and  wetting  of  the  polymer 
surfaces.  Once  good  contact  has  been  achieved  the  effect  of  pressures  from  0.3-2.8  MPa 
(50  to  400  psi)  appear  negligible.  At  very  high  pressures,  the  glass  transition  temperature 
of  PS  increases  by  about  20°C  per  kilobar  and  this  could  have  a  substantial  effect  on 
welding  rates  at  low  weld  temperatures. 

Stacer  and  Schreuder-Stacer  examined  the  effect  of  contact  pressure  on  the  autohesion 
of  PIB  (95).  Fig.  23  shows  their  data  for  PIB  surfaces  contacted  at  progressively  larger 
pressures.  The  data  were  obtained  as  described  previously  for  Fig.  16  and  the  curves  were 
shifted  using  the  same  WLF  coefficients.  Their  data  shows  that  as  the  pressure  increases, 
Ga  data  appears  to  shift  to  shorter  times  and  the  welding  time  too,  decreases  with  no 
discernible  effect  on  equilibrium  Goo  values.  They  suggest  that  the  data  support  the 
contact  mechanism  for  bond  formation  since  the  diffusion  process  would  be  independent  of 
pressure  in  this  pressure  range. 
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pig.  22.  Gjc  vs.  contact  pressure  for 
welding  of  polystyrene  at  constant  time 
of  30  min  and  temperature  of  1 15°C. 
(McGarel  and  Wool). 
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Fig.  23.  Effect  of  contact 
contact  pressure  on  the  autohesion 
curve  of  PIB.  (Stacer  and 
Schreuder-Stacer). 


Fatigue  of  Welded  Interfaces 

The  effect  of  welding  conditions  on  the  fatigue  crack  propagation  (FCP)  rate,  da/dN, 
along  the  welded  PS/PS  interface  was  investigated  by  Yuan  and  Wool  (98).  These  experi¬ 
ments  provide  information  on  the  expected  lifetime  of  welded  or  laminated  parts  as  a  func¬ 
tion  of  the  time  of  welding.  The  basic  question  being  addressed  is  to  determine  how  the 
FCP  rate  of  a  partially  healed  interface  compares  with  that  of  the  fully  healed  state. 

Double  cantilever  beam  (DCB)  specimens  as  shown  in  Fig.  24,  were  machined  from  care¬ 
fully  welded  plates  of  polystyrene.  A  thin  Al-foil  was  placed  in  the  interface  at  one  end  of 
the  DCB  to  promote  the  initial  crack  start.  Aluminum  blocks  with  pin  grips  were  glued  to 
the  DCB  ends  (Fig.  24).  The  fatigue  tests  were  conducted  in  displacement  <5,  control  mode 
with  a  haversine  wave  form  at  frequency  1  Hz  and  ratio  R  =  ^min/^max  =  0.33. 
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Fig.  24.  DCB  specimen  for  the  fatigue  testing  of  welded  polystyrene  plates  with  b  =  12.7 
mm,  d  =  12.7  mm  and  2h  =  6.6  mm.  (Yuan  and  Wool). 

The  crack  length  a,  was  measured  by  a  travelling  microscope  and  the  applied  strain 
energy  release  value  Gj,  was  calculated  using  the  following  expression  by  Wang  (99), 

G}  =  9(El/a4b)#[l  +  15/4a(<5/a)  (90) 
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+  (189/64  a2  +  6/7)(<3/a)2] 

where  E  is  the  modulus,  I  is  the  moment  of  inertia  and  a  =  (d  +  h/2)/a.  The  FCP  rate  is 
usually  a  power  law  in  aG  and  can  be  described  by 


da/dN  =  A4Gm  (91) 

where  A  and  m  are  constants.  In  this  expeiiment,  the  da/dN  points  were  taken  at  the 
proper  crack  length  to  ensure  applicability  of  beam  theory  and  linear  elastic  fracture 
mechanics 

The  FCP  rate  was  determined  at  each  welding  time  and  was  evaluated  at  a  constant 
Gt  =  240  J/m2.  The  results  are  shown  in  Fig.  25  for  welding  temperatures  of  108°C  and 
T  w  C  using  PS  samples  with  Mn  =  133,000  and  Mw  =■  303,000.  A  strong  effect  of  welding 
time  was  noted  and  the  results  can  be  expressed  as  a  power  law, 


(da/dN)a0  =  A(T)t-13  (92) 

where  A(T)  is  a  temperature  dependent  prefactor.  Thus,  as  the  strength  of  the  interface 
increases  with  welding  time,  the  FCP  rate  decreases  with  a  strong  power  law  dependence. 
The  exponent  of -1.3  compares  with  a  theoretical  value  of -1.25  predicted  using  a  chain 
disentanglement  model  (40,  50)  such  that 

(da/dN)AQ  a  (93) 

This  relation  suggests  that  at  t  =  Tr,  or  when  the  virgin  state  is  obtained  that  the  molecular 
weight  dependence  of  the  FCP  is  determined  by  (50), 

(da/dN)AG  a  M'5/2  (94) 

Fig.  26  shows  da/dN  vs.  molecular  weight  data  of  Hertzberg,  et  al.  (100)  and  the  strong 
molecular  weight  dependence  is  similar  to  that  predicted  by  Eq.  94.  It  is  interesting  to  note 
that  when  Berger  and  Kramer  (51)  examinecLtlie  craze  thickening  velocity  vc,  as  a  function 
of  molecular  weight,  they  found  that  vc  a  M'-V  .  The  agreement  between  these  results 
may  be  reasonable  under  stable  crack-craze  growth  conditions. 

The  effects  of  a  partially  welded  interface  on  the  FCP  rate  are  cons.derable  as  shown  in 
Fig.  25.  Many  plastic  parts  are  designed  with  regard  to  M*  such  that  Gjc  (see  Fig.  21)  is  a 
maximum  at  minimum  melt  viscosity.  However,  it  is  clear  that  enhanced  fatigue  lifetime 
can  be  obtained  with  M  >  M*,  as  can  be  deduced  from  Fig.  26.  These  results  are  of 
importance  for  composite  lamination  of  thermoplastic  matrices.  We  propose  further 
studies  of  fatigue-healing  problems. 

Fatigue  crack  healing  is  also  an  important  concept  for  fiber  filled  composites. 
Klosterman  and  Wool  (101)  found  that  damage  involving  fiber-matrix  failure  and  matrix 
cracking  in  glass  fiber  filled  thermoplastics  could  be  healed  at  temperatures  above  Tp.  As 
a  result,  the  FCP  rate  decreased  with  increasing  healing  time  in  the  damage  zone  unfil  the 
virgin  state  FCP  rate  was  obtained.  It  was  also  found  that  complete  crack  healing  could  be 
obtained  with  highly  cross-linked  epoxies.  These  studies  suggest  that  the  lifetime  of 
complex  composites  in  practice  could  be  considerably  increased  by  periodic  healing 
treatments.  The  possibility  also  exists  that  specific  microwave  and  other  treatments  could 
be  used  to  excite  healing  processes  in  specific  regions  of  the  composite. 


dA/dN  (  mm/cycle 
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Fig.  25.  Fatigue  crack  propagation 
rate  of  polystyrene  vs.  welding  time  at 
constant  G  =  240  J/m2.  Mn  =  133,000 
and  Mw/Mn  =  2.3.  (Yuan  and  Wool). 


Fig.  26.  Fatigue  crack  propagation  rates 
vs  molecular  weight  of  PVC  (Rimnac,  et  al.). 
The  solid  line  is  given  by  disentanglement 
theory  by  daidn  =  alt Ml  ' z  -  Mcp  '■ z)  where 
a  =  4.5  x  10°  mm  z  cycle  and  Mc,-  = 
60,000  is  the  molecular  weight  at  which 
critical  crack  propagation  occurs  at  K  =  0.7 
MPa  Vm. 


Incompatible  Amorphous  Interfaces 

With  immiscible,  incompatible,  amorphous  polymers  the  equilibrium  structure  of  the 
interface  is  determined  thermodynamically  by  the  balance  of  (1)  the  decrease  of  free 
energy  due  to  entropy  relaxation  of  the  surface  chain  configurations  and  (2)  the  positive 
enthalpy  of  mixing  produced  by  limited  diffusion  to  a  distance  d».  Since  do,  ~  X'1/" 
according  to  Helfand,  oi  al.  (24,  86)  and  the  Flory-Huggins  interaction  parameter,  X. 
behaves  with  temperature  approximately  as  X  «  1/T,  we  expect  for  fracture  dominated  by 
chain  pullout  that  Gjr  ~  loo  ~  doo2  and  that  the  equilibrium  strength  behaves  with  welding 
temperature  as,  (22,  77) 


Goo  -  T.  (95) 

The  equilibrium  fracture  energy  for  several  incompatible  interfaces  was  obtained  by 
Foster,  Willett,  et  al.  using  wedge  cleavage  methods  (77,  87. 102).  These  included 
PS/PMMA,  PSAN/PMMA,  PC/PSAN.  Gjr  was  found  to  increas'  with  temperature  in 
reasonable  agreement  with  Eq.  95  for  PSAN/PMMA  and  PC/PSAN  interfaces  but  was 
fairly  independent  of  temperature  for  PS/PMMA  up  to  170  C..  Typical  measured  fracture 
energies  of  order  40-200  J/m2  were  greater  than  either  the  initial  wetted  energy  (about  5 
J/m2)  or  the  average  surface  energy,  T  <  1.0  J/m2,  indicating  that  limited  difmsion  had 
taken  place.  The  dynamics  of  welding  incompatible,  amorphous  interfaces  was  found  to  be 
complex  and  will  not  be  presented. 
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The  equilibrium  strength  of  incompatible  interfaces  occurs  with  full  contact  of  the 
surfaces  and  appears  to  be  dependent  on  the  extent  of  the  limited  interdiffusion  which  is  a 
function  of  temperature,  molecular  weight  and  chemical  compatibility.  For  example,  in 
Fig.  14,  the  SBR/Butyl  interface  in  the  autohesion  experiment  is  incompatible  and  low 
strength  is  achieved  compared  to  the  symmetric  SBR/SBR  and  Butyl/Butyl  interfaces  even 
though  good  contact  has  been  achieved  at  very  long  times.  This  observation  strongly 
supports  the  diffusion  mechanism  of  strength  development  at  amorphous  polymer 
interfaces. 

The  magnitude  of  the  equilibrium  interpenetration  distance,  d°°,  can  be  approximated 
from  the  mechanical  data  using  the  following  scaling  relations, 


[doo/SAP  =  Gk(A/B)/GIc(A/A) 
or, 


(96) 


[<WSA]  -  o(A/ B)/ct(A/A)  (97) 

where  SA  is  the  radius  of  gyration  of  the  pure  A-component,  Gjc(A/B)  and  a(A/B)  are  the 
equilibrium  fracture  properties  of  the  incompatible  pair  and  Gjc(A/A)  and  o( A/A)  are  the 
virgin  state  fracture  properties  of  the  A  component.  It  is  assumed  that  the  pure  polymers  A 
and  B  have  the  same  virgin  state  mechanical  properties  and  molecular  weights.  For 
example,  the  butyl  elastomer  in  Fig.  14  has  a  viscosity  average  molecular  weight,  Mv  = 
225,000  such  that  S  «  200  A.  If  a(A/B)/a(A/A)  «  100/400,  then  from  Eq.  97,  we  obtain 
doo  ~  50  A. 

For  PS/PMMA  interfaces,  Foster  and  Wool  (87),  using  wedge  cleavage  methods 
obtained  constant  values  of  GIc  =  45  J/m2  for  welding  at  temperatures  between  100  and 
170°C.  If  we  assume  that  for  polystyrene  with  M  =  250,000  and  S  =  138  A  that  GIc(A/A) 

=  1000  J/m2,  as  shown  in  Fig.  21,  then  from  Eq.  96,  d»  =  29  A.  The  Flory-Huggins 
interaction  parameter  X,  is  related  to  d®  by  X  =  0.66  (b/d°°)2,  where  b  ~  6  A  is  the 
statistical  segment  length.  If  d°°  =  29  A,  then  X  «  0.028  for  the  PS/PMMA  interface.  The 
temperature  independence  of  the  fracture  energy  for  PS/PMMA  interfaces  are  consistent 
with  the  very  weak  temperature  dependence  of  X  measured  by  Russell,  et  al.  using  neutron 
scattering  techniques  (103). 


Incompatible  Crystallizable  Interfaces 

Strength  development  at  incompatible,  crystallizable  interfaces  such  as  polyethylene 
with  polypropylene,  PE/PP,  can  occur  by  different  mechanisms  than  those  discussed  above. 
In  the  melt,  limited  diffusion  to  a  depth,  do,,  first  occurs.  When  the  interface  cools,  crystal¬ 
lization  occurs  first  in  the  PP  side  due  to  the  difference  in  melting  points.  Random  nuclea- 
tion  of  spherulites  near  the  interface  coupled  with  the  volume  contraction  due  to  crystal¬ 
lization  breaks  up  the  interface.  PE  melt  riverlets  are  drawn  into  the  PP  side  which 
subsequently  crystallize  and  effectively  mechanically  interlock  the  interface.  This  mech¬ 
anism  has  been  discussed  by  Galeski  (88,  89)  in  terms  of  "local  crystallization"  creating 
"influxes"  across  the  interface  and  has  been  investigated  by  Yuan  and  Wool  (77,  90). 

Fig.  27  shows  a  PE/PP  interface  during  fracture  by  wedge  cleavage  methods.  The 
interface  was  crystallized  in  a  manner  to  optimize  the  mechanical  interlocking  effect.  The 
fibrillar  structures  in  the  fractured  interface  are  due  to  plastic  deformation  of  the  PE 
influxes.  The  fracture  stress  was  measured  as  a  ~  10  MPa  which  is  more  than  half  the 
virgin  strength  of  PE.  In  the  absence  of  influxes,  the  interface  is  extremely  brittle  with  a  < 
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1  MPa.  Fig.  28  shows  a  magnified  view  of  the  PP  side  obtained  by  SEM  (90).  The  PE 
fibrillar  remnants  are  clearly  seen  in  the  interstices  of  the  PP  spherulites.  In  addition  to  the 
influxes,  PP  spherulites  which  nucleate  on  the  interface  plane  entrap  PE  melt  on  a  finer 
scale  between  the  crystalline  lamellar  structures  as  seen  by  the  "white  decoration"  of 
plastically  deformed  PE  on  several  PP  spherulites. 


Fig.  27.  Fracture  of  a  PP/LLDPE  incompatible  semi-crystalline  interface.  The  fibrils  on 
the  PE  top  surface  are  due  to  pullout  of  PE  material  which  was  incorporated  into  the  PP 
side  during  crystallization  (Yuan  and  Wool). 


Fig.  28.  (a)  Electron  micrograph  of  a  PE  fracture  surface  from  a  PP/LDPE  interface. 
Plastically  deformed  PE  fibrils  are  observed  in  addition  to  circular  imprints  of  PP 
spherulites  on  the  PE  side,  (b)  PP  fracture  surface  from  a  PP/LDPE  interface  shows  a 
fractured  PE  fibril  remnant  and  PP  spherulites  (Yuan  and  Wool). 

Thus,  for  incompatible,  crystallizable  interface,  while  melt  diffusion  produces  an 
equilibrium  interface  of  thickness  50  A,  the  localized  crystallization  mechanism  produces 
influxes  to  depths  of  several  microns  which  impart  significant  strength  to  the  interface. 
Delicate  control  of  the  processing  history  is  necessary  to  optimize  this  effect  as  it  is  easy  to 
make  a  weak  interface  with  incompatible  polymers.  Methods  of  optimizing  the  strength 
have  been  discussed  (90).  The  same  mechanism  also  exists  at  the  interface  between 
symmetric  crystallizable  interfaces  such  as  PE/PE  welding  but  plays  a  secondary  role  to 
diffusion  and  homogeneous  crystallization. 
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Lap  Shear  Welding 

A  lap  shear  method  as  shown  in  Fig.  29  was  used  to  study  strength  development  during 
welHina  0f  PS  interfaces  by  Kline  and  Wool  (91).  The  surfaces  had  previously  not  been  in 
contact  and  care  was  taken  to  ensure  rapid  wetting  of  the  interface.  Lap  shear  samples 
were  welded  for  a  given  time  at  constant  temperatures  above  To  and  fractured  at  room 
temperature.  Fig.  30  shows  a  plot  of  the  (average  of  15  samples)  shear  failure  stress,  rf,  vs 
time  to  the  fourth  power  at  each  constant  temperature.  The  linear  response  of  rf  with  t1' 4 
is  seen  at  each  temperature  and  provides  further  support  for  diffusion  controlled  strength 
at  polymer-polymer  interfaces. 


Fig.  29.  Lap  shear  joint  method  used  to  test  weld  strength  of  polymers  (Kline  and  Wool). 

The  change  in  slope  with  temperature  was  found  to  be  describable  by  the  William- 
Landel-Ferry  (92)  theory  of  thermal  activation  near  Tg.  At  T  =  113.5°C,  the  activation 
energy,  Ea  =  93.2  kcal/mole  is  predicted  from  a  WLr  analysis  using  the  parameters  for  PS, 
ci  =  13.7,  C2  =  50  and  Tg  =  100°C,  which  compares  with  Ea  =  96.1  kcal/mole  obtained 
from  a  plot  of  the  superposition  shift  factor  vs.  1/T  in  Fig.  30. 


The  time  to  achieve  complete  healing,  T^,  can  be  estimated  from  the  upper  plateau  in 
Fig.  30.  At  391.3°K  (118.1°C),  Tj-,  *  256  min  for  the  polydisperse  PS  samples  with  = 
142,000  and  Mw  =  262,000.  This  time  can  be  compared  with  related  times  measured!  by 
other  techniques.  Lee  and  Wool,  (14)  using  FUR  dichroism  techniques,  measured  the 
relaxation  time,  of  uniaxially  oriented  monodisperse  PS  with  Mn  =  233,000  as  Tr  =  215 
min.  Kramer,  et  al.  (93)xeport  a  self-diffusion  coefficient  for  monodisperse  PS  with  Mn  = 
255,000  as  D  =  5.8  x  10' 1 '  cm2/sec  at  125°C.  Shifting  to  118°C  using  the  relation,  (93) 


log  D/T  =  A  -  B/(T-Too), 


(98) 
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Fig.  30.  Shear  stress  vs.  for  lap  welding  of  polystyrene  at  several  temperatures  (Kline 
and  Wool). 

where  B  =  710,  Too  =  49°C  and  A  =  -9.49,  then  Diig°c  =  6.42  x  10'^  cm2/sec.  Using 
Eq.  76,  the  time  required  for  the  center-of-mass  to  diffuse  a  distance  equal  to  the  end-to- 
end  vector  is  obtained  as  Tr  =  1000  min  which  is  about  4  times  longer  than  the  time 
required  to  achieve  complete  strength  of  the  lap  shear  joint  with  Mw  =  262,000.  On  the 
other  hand,  using  Tr  ~  M3  and  D  =  M-2,  the  relaxation  time  for  the  chains  with  Mn  = 
143,000  at  118°C  is  about  180  min  which  is  in  close  agreement  with  the  lap  shear  result. 
Because  of  the  ramified  nature  of  polymer  interfaces  with  holes  and  statistical  irregularities 
as  shown  in  Fig.  7,  the  mechanical  healing  times  are  expected  to  be  a  little  longer  (but  of 
the  same  order  of  magnitude)  than  the  relaxation  times  determined  from  molecular 
dynamics  and  diffusion  studies.  It  should  be  noted  that  when  evaluating  healing  times  from 
the  reptation  theory  that  at  t  =  Tf,  the  center-of-mass  of  chains  in  the  bulk  diffuses  a 
distance  related  to  the  end-to-end  vector  (Eq.  76),  while  the  average  monomer  interpene¬ 
tration  distance  at  the  polymer  interface  is  equivalent  to  the  radius  of  gyration  (Eq.  17). 

POLYMER  PROCESSING  WELD  LINES 


Compression  Molding 

As  a  final  example  of  welding  we  consider  the  compression  molding  of  resin  pellets 
where  the  interfaces  and  weld  lines  form  when  the  pellet  surfaces  contact  in  the  mold  as 
shown  in  Fig.  31.  With  increasing  time  and  pressure,  the  interfaces  wet,  diffusion  occurs 
and  the  strength  of  the  sample  increases  to  its  virgin  strength  provided  that  sufficient  time 
for  diffusion  is  allowed. 


Fig.  31.  Schematic  of  powder/pellet  interface  welding  during  polymer  melt  processing. 
With  time,  temperature  and  pressure,  the  pellet  interfaces  heal  to  form  the  plastic  product. 

Wool  and  O'Connor  (1)  studied  the  time-dependence  of  compression  molding  using  PS 
pellets  of  the  same  molecular  weight  as  described  above  in  the  lap  shear  experiments.  PS 
plates,  152  by  152  mm,  were  compression  molded  in  a  Carver  Press  using  a  window  frame 
3.18  mm  thick.  Each  plate  was  molded  at  127°C±  6°C,  with  pressures  of  5.5  MPa  (800  psi) 
for  varying  processing  times,  tp.  The  plate  was  withdrawn  from  the  mold,  quenched,  cut 
into  tensile  "dogbone"  sampled  and  fractured  at  room  temperature  in  uniaxial  tension. 
Weak  interfaces  within  the  sample  provided  a  source  of  crack  initiation.  Fig.  32  shows  the 
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fracture  stress  (average  of.6,  samples)  plotted  vs.  the  fourth  power  of  the  processing  time. 
The  prediction  of  a(t)  ~  t1/4  is  again  noted  to  be  satisfactory  despite  the  crudeness  of  the 
experiment.  The  processing  time  to  ieach  the  virgin  state  near  ct0 0  ~  41.4  MPa  (6,000  psi) 
was  between  10  and  20  min  judging  by  the  saturation  of  the  strength  at  longer  times.  Using 
Kramer's  diffusion  data  at  125°C,  the  reptation  time  for  pellets  with  Mn  =  143,000  is  about 
20  min  which  is  in  reasonable  agreement  with  the  observed  processing  time. 


Fig.  32.  Fracture  stress  vs.  processing  time  (t1/4)  for  compression  molding  of  polystyrene 
pellets.  Dogbone  samples  were  prepared  from  compression  molded  plates  (Wool  and 
O'Connor). 

The  molecular  weight  dependence  of  the  virgin  strength,  CToo,  of  compression  molded 
PS  pellets  has  been  determined  by  McCormick,  et  al.  j$4)  using  similar  sample  preparation 
and  fracture  conditions.  Fig.  33  (1)  shows  a vs.  for  monodisperse  molecular  weight 
PS  using  data  from  reference  94.  The  data  are  described  by 

Ooo  =  0.23(M1/2  -  Mi1/2)  MPa,  (99) 

where  the  intercept  at  a  =  0  is  VMj  =  245.  The  strength  is  seen  to  increase  linearly  in 
accordance  with  virgin  state  predictions  up  to  M  «  160,000  after  which  it  decreases  slightly 
with  increasing  M.  The  data  was  obtained  from  samples  compression  molded  at  a  constant 
time  of  1  min  at  170°C  for  all  materials.  We  have  suggested  that  under  the  processing 
conditions  used  that  samples  with  M  <  160,000  had  reached  their  virgin  state  at  t  =  1  min 
but  that  samples  with  M  >  160,000  were  not  completely  healed  at  1  min  and  the  fracture 
stress  should  decrease  as  observed  according  to  a  =  M'1/4.  The  reptation  time  for  M  = 
160,000  is  about  1  sec  at  170°C  and  the  time  of  1  min  used  in  this  experiment  presumably 
includes  warm  up  time  in  the  mold  in  addition  to  diffusion  times.  However,  the  maximum 
in  strength  with  increasing  M  at  constant  contact  time  is  expected  and  is  similar  to  the  tack 
and  green  strength  data  presented  in  Fig.  15. 

Strength  development  at  internal  weldlines  in  injection  molded  specimens  are  expected 
to  behave  similarly  although  some  complications  are  expected  with  interfaces  involving 
oriented  molecules  at  the  weldline.  However,  it  is  interesting  to  note  that  orientation 
relaxation  occurs  at  about  the  same  rate  as  normal  diffusion  (14).  Orientation  relaxation 
may  contribute  to  interdiffusion  at  the  weldline  in  a  similar  manner  as  for  unoriented 
surfaces  but  problems  are  anticipated  with  "barrier"  effects  produced  by  oriented  layers 
which  could  result  in  anisotropic  diffusion  rates. 


Fig.  33.  Tensile  strength  vs.  square  root  molecular  weight  of  compression  molded 
polystyrene  (pellets).  Molding  occurred  at  170°C  for  1  minute.  Monodisperse  molecular 
weights  pellets  were  used.  Data  of  McCormick,  et  al. 


Relation  to  Vibration  Welding 

Vibrational  welds  can  be  affected  in  thermoplastics  over  a  frequency  range  25-400  Hz, 
as  reviewed  by  Stokes  (104-107).  Typical  welding  pressures  are  130-2000  psi  (1-14  MPa) 
with  weld  amplitudes  of  0.3  to  1.6  mm.  The  welding  process  can  be  divided  into  four 
phases  which  are  independent  of  the  thermoplastic  as  follows.  In  Phase  1,  the  weld 
interface  is  heated  by  Coulombic  friction  which  ends  when  the  interface  begins  to  melt. 

The  penetration,  which  is  the  moving  together  of  the  separate  pieces  due  to  lateral  flow  out 
of  the  interface,  is  essentially  zero  during  this  phase.  With  continued  energy  input,  the 
interface  begins  to  melt  and  flow  and  penetration  commences  in  Phase  2  up  to  a  steady 
state  lateral  flow.  Phase  3  represents  the  steady  state  flow  condition  during  which  the 
penetration  increases  linearly  with  time  such  that  the  material  is  melting  at  the  same  rate 
as  it  is  flowing  out  of  the  lateral  surfaces.  The  molten  film  is  of  the  order  of  0.1  to  0.23  mm 
and  the  temperature  is  very  close  to  ambient  within  about  1.0  mm  of  the  solid-liquid 
interface  (105).  Phase  4  occurs  when  the  vibratory  motion  is  stopped,  and  the  interface 
solidifies. 

The  role  of  diffusion  in  vibrational  welding  is  apparent  in  Phases  2-4.  However,  since 
the  diffusion  necessary  to  promote  maximum  strength  involves  distances  of  the  order  of  100 
A,  the  large  amount  of  local  melting  (0.1  mm)  and  flow  may  be  unnecessary.  Obviously, 
the  extent  to  which  Phase  3  steady  state  flow  is  sustained  determines  the  extent  of  lateral 
flow  and  local  deformation  of  the  weldline.  For  each  thermoplastic,  the  temperature  rise 
at  the  end  of  Phase  2  could  be  used  to  determine  the  healing  time  from  the  diffusion 
analysis  and  this  could  be  used  to  minimize  Phase  3  flow  processes.  For  optimal 
component  integrity,  we  would  like  to  minimize  contributions  from  Phases  2  and  3  by 
heating  the  interface  in  Phase  1  to  cause  melting  and  solidify  in  Phase  4.  The  pressure 
history  could  also  be  varied  to  promote  melting  in  Phase  1  but  with  reduced  pressure  in 
Phases  2-4  to  minimize  flow.  The  reduced  pressures  used  in  vibrational  welding  would  not 
affect  the  interdiffusion  process  once  good  contact  had  been  achieved.  The  fact  that  flow 
occurs  during  Phase  3  at  the  interface  presents  new  opportunities  for  mechanical 
interlocking  and  this  could  be  important  for  incompatible  interfaces. 


CONCLUDING  REMARKS 


In  this  Report  on  welding  mechanisms,  we  have  examined  strength  development  in 
terms  of  1)  the  structure  of  the  welded  interface,  2)  the  microscopic  deformation 
mechanisms  and  3)  the  fracture  mechanics  of  the  weld.  Attempts  were  made  to  relate  the 
structure  to  strength  using  entanglement  concepts  and  the  dynamics  of  random-coil  chains. 
The  status  of  the  description  of  the  symmetric  interface  is  such  that  we  have  significantly 
detailed  information  from  the  concentration  profiles  as  shown  in  Table  1,  provided  the 
reptation  model  is  the  correct  dynamics  model.  It  is  likely  that  more  detailed  models  will 
evolve  from  current  research.  However,  if  one  had  to  hypothesize  at  the  nature  of  new 
models,  it  is  most  probable,  based  on  the  entanglement  constraints  to  motion  in  polymer 
melts,  that  the  framework  of  the  reptation  model  would  still  be  visible.  The  fractal  nature 
of  the  diffusion  front  is  generic  to  diffusion  processes  and  provides  useful  physical  insight 
into  the  structure  and  properties  of  interfaces. 

The  relation  between  structure  and  strength  is  complex  and  necessitates  the  use  of  new 
entanglement  concepts  in  microscopic  deformation  models.  The  latter  models  and  several 
simplifying  assumptions  used  in  this  Report  in  relating  structure  to  strength  may  be  a  cause 
of  concern  to  the  reader.  This  self-criticism  points  to  the  gap  in  our  understanding  of 
entanglement  concepts  and  mechanical  properties  at  the  microscopic  level.  For  example, 
the  3.4  exponent  for  the  molecular  weight  dependence  of  the  zero-shear  viscosity  remains 
an  unsolved  problem  and  impacts  on  our  confidence  to  relate  molecular  dynamics  to 
mechanical  properties.  It  was  the  hope  of  this  Report  to  present  an  analysis  of  strength 
development  at  polymer-polymer  interfaces  in  a  manner  to  indicate  both  their  interesting 
complexity  and  suggest  pathways  for  solutions  to  interface  problems. 
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